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TRAKSLATOR'S  PREFACE 

As  research  in  dynaiillcal  weather  prediction  proceeds,  increasing 
attention  is  tile Ing  given  to  the  mathematical  pi>3blemB  encountered  in 

i;  ,  ,y  ■  ; 

the  computation  of  humarical  aolutiona  oi  the  dynamical  equatlonB  for 
the  atmosphexe;  the  meteorologist's  attention  Is  thus  drawn  to  problems 
of  numerical  ahelysls .  Basic  to  much  of  this  work  are  the  classical 
existence  and  unlqiieness  proofs  of  Oourant,  Friedrichs  end  Levy  for 
elliptic  and  hyperbolic  partial  difference  equations.  The  present 
translation  has  been  prepared  In  order  to  fill  the  heed  for  an  English 
version  of  this  paper,  |ud  thereby  to  make  this  fundamental  work  luore 

"  II  '  ' 

readily  available  and  useful  to  the  meteorological  profession  In 
paxtlcu^^n  After  this  translation  was  completed,  however, ,,mjf  atten¬ 
tion  was  dravm  to  a  previous  English  translation  prepared  by  M.  H.  Rand 
in  19^6  at  „t^e  British  Atonic  Energy  Besearch  Establishment,  Harwell. 
Tha  ijtoitqwrlBOD  is  in  general  reassuring,  and  a  number  of  errors  noted 
in  Band's  translation  have  here  been  avoided. 

partl^ulfr  .oasv  has  been  taken  in  the  present  translation  to 
preseive  the  spirit  of  the  original  Oeinian  as  closely  as  possible.  In 
this  connection  Miss  OberlEnder,  of  .the  Departawnt  of  Meteorology, 
UVC.L.A.,  contributed  significantly,  as  well  as  assisting  In  the 
preparation  of  the  final  manuscript.  The  figures  have  been  redrafted 
and  enlarged  to  improve  their  readability,  although  no  other  changes 
c^f  format  have  been  made. 

M.  Lawrence  Gates 
Project  Director 


ON  THE  PARTIAL  DIFFERENCE  EQUATIONS  OF  MATHEMATTOAT.  PRYHICB 

By 

R.  Courant,  K.  FTledrichs,  and  H.  Iicwy 

In  dfsnling  with  classical  linear  differential  egjiations,  if 
one  replaces  the  differential  quotients  by  difference  quotients 
defined  over  some  — •  let  us  assume  rectangular  —  grid,  one  obtains 
algebraic  problems  of  a  much  more  transparent  structure,  Th^ 
following  paper  wi]:l  undertake  an  elementary  discussion  of  these 
algebraic  problems  isnd  'ill  principally  discuss  the  behavior  of  the 
solution  as  the  mesh  width  of  the  grid  goes  to  zero.  To  these  ends 
we  limit  ourselves  mainly  to  the  simplest,  but  nevertheless  typical, 
cases  which  we  handle  in  such  a  manner  that  the  applicability  of  our 
methods  to  more  general  difference  aquations  with  arbitrarily  many 
Independent  variables  will  be  clear. 

Corresponding  to  the  familiar  problems  of  differential  aqua¬ 
tions,  we  will  ti-eat  boundary  value  and  eigenvalue  problems  for 

elliptic  difference  equations,  and  the  Initial  value  problem  for 

ji  '  "  ■ 

hyperbolic  or  parabolic  difference  equations.  We  shall  prove  with 

sloqtle  typical  examples  that  the  limiting  process  Is  always  possible 
l.e.,  that  the  solutions  of  the  difference  equations  converge  toward 
the  solutions  of  the  corresponding  dlffei'ential  equations;  indeed  we 
shall  find  that  for  elliptic  equations  In  general,  the  difference 
quotients  of  arbitrarily  high  order  tend  toward  the  corresponding 
differential  quotients.  The  existence  of  a  solution  to  the  differ¬ 
ential  aquation  is  nowhere  assumed  and  indeed  we  obtain  through 


But  while  for  the 


the  limiting  process  a  simple  existence  r.roof.-'' 
elliptic  case  convergence  is  the  rule  Independently  of  the  choice 
of  grid,  ve  will  show  that  in  the  case  of  the  initial  value  problem 
for  hyperbolic  equations  convergence  can  be  assumed  only  when  the 
ratio  of  ffrid  mesh  sizes  in  different  directions  satisfies  certain 
inequalities,  which  ip  turn  are  determined  by  the  ptJbiLiou  of  the 

'r/ 

cbaracterieticui  relative  to  the  grid. 

We  take  as  a  typical  example  of  the  elliptic  case  the 
boundary  value  problem  of  potential  theory.  The  dependence  of  its 
solution  on  the  solution  of  the  correopondlng  difference  equation 
has  been,  moreover,  extensively  treated  during  the  i>aBt  few  ^ars.^ 
Of  course,  in  our  case,  in  contrast  to  the  previous  work,  the 
particular  special  properblea  of  the  potential  equation  will  not 
be  explicitly  used  so  that  the  application  of  our  method  to  other 
problems  is  not  overlooked. 


1  Our  method  Of  proof  may  ba  extended  without  diffloull^  to 
cover  the  boundary  and  eigenvalue  problema  for  arbitrary  linear 
elliptic  dlffei*ential  equations,  and  the  initial  value  problem  for 
arbitiary  i.inear  hyperbolic  differential  eqpetiona. 

,1,  le  Roux,  SuT  le  probliae  de  DlrlcUlet,  Joum.  de  aathem. 
pur,  at  appl . .  (6)  10  1^191^)1  p.  I89,  R.  G.  D.  Richardson,  A  New 
MetHojT'iTi  Boi,indary  frbbl^pBS  for  Differential  Equations,  Trees.  Amer. 
Math,  80c.  16  (1917);  ,P‘f  H.  B,  Phillips  and  N.  Wiener,  Nets 

«jd  the  Dlrioklet  Problem,  publ,  of  th«i  Maas.  Institute  of  Technology 
(1985),  ( Translator 'b  note?  see  J.  Math,  Phyo..  fi;105-12U  (1985)). 

Unfortunately  these  papers  were  mlaaed  by  the  first  of  the 
present  thi'se  autluirs  (Courant)  in  the  writing  of  his  note  "Zur 
Theorio  d«r  partiallen  Dlifarenzengleidtungen",  a(?tt.  Nachr, .  25>  X. 
1Q?5,  which  the  prfiBcnt  Work  extendn. 

Compare  furt'aftr;  L.  Lusternlk,  "(jbor  etnige  Anvendungen  der 
direkten  Methcdan  in  der  Variatlonsrechnung",  loculi  de  1^  Societe 
Mathsm.  rle  Woncou,  I926,  G.  Bouligand,  ”3ur  Ie"proFlcme  do  Dirichlet" 
/cmT  de  la  bqc  .  polon.  dg  laathem. ,  4  Krakau,  1986, 

For  the  of  difference  exqji’eeBlons ,  and  for  related  works, 
see  R,  Courant,  "Uber  dlrekte  Mefchodeu  in  der  Varlatlonsrechnung", 
Math.  A.iin-alan ,  97,  p,  'fll  and  the  literatr^ire  cited  there. 
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In  addition  ■fco  the  stated  priiAGipal  part  of  the  work,  ws  shall 
append  an  elementary  algebraio  dischsaion  of  the  boundary  value  problem 
for  the  elliptic  equations  connected  with  the  well-known  problem  of 
random  walks  arising  in  statistics. 


1.  fHk  kLLlpTiiJ  CASE 
$  1.  Preliminary  Remarks 

1.  Definitions .  In  the  plane  with  rectangular  coordinates 
X,  M  we  consider  first  of  all  a  sqjiare  grid  of  pointa;  of  mesh  width 
(t.>0  ,  all  points  having  the  ooordinateB“  X  -  ,  ysfnk  , 

h, (n  a  «>,  ±1,  ta,  . 

Nov  let  G  be  a  region  of  the  plane  bounded  by  a  continuous 
closed  curve  that  is  free  of  double  points ,  Then  the  corresponding 
grid  region  -  vhloh  Is  uniquely  determined  for  eufflolently  small 
mesh  width  -  sh^ii  consist  of  all  those  grid  points  that  lie' in  G  end 
which  can  be  Joined  to  a  fixed  or  preseribed  grid  point  of  G  %  & 
connected  series  of  grid  points.  We  denote  as  a  eonneoted  series  of 


grid  points  a  sequence  of  poihtii  auoB  tha^  WMb  point  in  Ipia ^ 
of  the  four  neighboring  pdints  of  the  folloving  point.  We  denote  ae 
a  boundary  point  of  G|^  a  point  vboee  four  neic^boring  points  do  not 


all  belong  to.  G|^ 


All  other  points  of  G|^  we  e*ll  Interior  points. 


We  shall  consider  functylons  U  $  Vt  •••  of  position  in  the  grid, 


i*e.,  functions  which  are  defined  only  for  grid  points.  We  shall  also 
denote  them  as  ( X  ,  y  ),  v(X  » y  )>  ■  •  •  •  their  forward  end 
backward  difference  quotients  we  employ  the  following  abbreviations: 


-k- 


S 


-^[K.(x,y)''M.(x-K.^=  u-  >  “j^u(x,y)“U(x,»j-K)j= 

Corresponcllngly  we  form  the  difference  quotienta  of  higher 
order,  eig., 

*  Wx5t  “  '*^3?x  *  ■^|'Wfx4K,y)“2-»tCx,^)  +  itCx-h,')>],  ctc 

2.  Dlfferenoe  Expreaslona  and  Green's  Trans formatlona.  We 
”got  the  slB^IeBt  general  view  or  linear  difference  expresaions  of  the 
second,  order  from  the  pattern  of  the  theory  of  partial  differential 
equations  if  we  form  from  two  ftmotions  w  and  V  and  their  forward 
dlfferenoe  quotients  a  bilinear  egression 

B(u,v)a  +  bw,fV,^+ 

tot u^v  +  f  i'kVyfi'MVy# 3tiv, 


where 


a s  aCKfi),  •  ^  ®  ^ 


are  functio&e  ia  the  grid. 


-  Pk)b  the  bilinear  e3®i»  of  the  first  order  we  derive  a 

■  -- «  '  //  : .  '  ,  . 

difference  eiqpreBBlpn  of  the  second  order  in  the  following  way:  ve 

form  the  sum 

4 


over  all  points  of  a  region  /;  in  the  grid  wherein  B(K,V’)  Is  set 

K 
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to  zero  for  t,he  dlfferenoe  (mottent  between  a  boundary  point  and  a 
point  that  does  not  belong  to  We  now  transform  the  sum  through 

partial  summation  (i.e.,  we  arrange  according  to  V ),  and  split  it 
Into  a  sum  over  the  set  of  interior  points  and,  a  aum  over  the 
set  of  boundary  points  .  Thus  we  obtain 

-  vL(u}-hS  VR(U).  (1) 

..  /k  ■; 

L(u)  is  the  linear  "diffei'ence  expression  of  second  order"  defljied 
for  all  interior  points  of  Gi^, 

LCuTs  C<mKH  +  (^'*dx)g  +  (cuy)y  +f«iuy)g 

-  --  4  (<“u)5(  •+  ' 

f((VC)  is  a  lineer  difference  expression  for  every  boundary  point 
whose  eicaet  form  we  shall  not  give  here. 

.  If  we  arrange  2.  SI  8('k>V')  with  respect  to  V- ,  we  get 

,  .  Iv  tiM/v)-k2lKS(ir).  (2) 

Gf,  Ofl  Hi 

M(V’)  ie  sailed  the  adjoint  difference  expression  of  L(  U)j  ve  have 

^  +  (l>v^)x  +  (^^^5 

;  +  (o«v)r+ ~  ‘  » 

where  $  ( v)  la  a  difference  expression  for  the  boundary  corresponding 
to  R(U). 

The  equations  (1),  (2)  and  the  following  equation 

h^51!E(vU'u)-'uM(v'))  4  k£I  (vR('u)-uS(v))»9  (3) 

K  , 

are  called  Green's  formulas.  " 
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The  Bimpleat  and  most  iraportont  case  occurs  if  the  hilinear 
form  iB  symmetric,  i.e.,  if  the  equations  b*  C  ,  0<=  )f  ,  ^ 

exist.  In  this  case  the  expression  L(K-)  la  Identical  with  its 
adjoint  M(V)j  we  therefore  call  it  the  self-adjoint  caee,  and  it  - 
is  readily  phtalnahle  from  the  quadratic  expression 

— -u  .  ' 

In  the  follovlng  ve  shall  limit  ourselves  mainly  to  expressions 
1^  (14,)  which  aws  self -adjbl^,  ^hftjoliaracter  of  the  diffaienee 
expression  depends  atovei  all  on  the  natural  of  those  term  of 

the  qpadratlo  form  vhleh  are  qjiadMtle  in  the  first  dif- 

feranoe  quotients.  We  call  this  part  of  B(tt,  l4)  the  eharacterist^ 


. :// 


p  (t4,  H  )  «  an  J 


form! 


Mcordii^y  as  f  (tiiH)  ia  (positive)  definite  or  indefinite  in  the 
differsnoe  qiixotlstits  ve  call  the  oorrsspqndlng  differenoe  expresilon 


or  l^ypsrt^  ■  . 

The-ditfsrsnaa  -sj^raasion'.;. 

With  Which  ve  shall  primarily  concern  ourselves  "in  tha  fellowlai 

‘i-i  ' 

paragraphs,  is  elliptic..  In  parfcloular,  it  is  obtained  from  the 


.u .  // - ^.^4 


quadratic  expression 


B  (U,'k  )  s  +  n  *  OP 


a  a 


I 

IL 
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asnce  the  corresponding  Green's  fonjiulas  are^ 


The  difference  expression  Att  «  obviously  the 

analogue  of  the  differential  ej^resalon  2>V^dX* for  a  fUactloh 
)  of  the  eontlnuoua  variables  X  and  u  .  tfrittea  out 
'  eompletely  the  difference  expression  reads 

Therefore  K  is  the  excess  of  the  arithmetic  mean  of  the 

,/  function  values  at  the  four  neighixring  polnte,  over  the  function 

.  •  ,  .  •  '  •<*"  ., 

■  value  at  .the, point  in' qjiestion,.-"  ■ 

'Gon^letely  sinaiar  Consideratiens  1^^^  to  linear  dlfferenoe 
:  r^A  ”  order  and  oorresponding  Green '^e  fonulas  if 

-  -y*  ^gin  witK  expressions  which  are A  fproed 

from  the  difference  ^otients  of  second  order.  We  content  ourselves 
with  the  example  of  the  difference  expression 

•»  ••  A.  '  '  '■ 

AAu  a  Sxsfx  ’^  '^^55  '  , 

^  The  boundary  expression  R(tL)  may  be  written  here  as  follows  t 
If  Ho,  H,  ,  be  the  function  values  at  the  boundary  point 

concerned  and  at  Its  .0  neighboring  points  ( i>  STS),  then 

R  (H)  »  ( 11,  -  ii>«a  )  . 


fi,'/  ■ 


■7" 


This  arises  from  the  quadratic  expression 


=:  ("U^x  +  ^yy)  == 


If  one  orders  the  sum 


li’  ZS 


with  respect  to  V ,  or  correspondingly  replaces  *14.  hy  the  expression 
AM,  in  equation  (5).  One  must  notice  however  that  in  the  expression 
AAu,  the  function  value  at  a  point  is  connected  with  the  function 
values  at  its  nelgh'boHng  points  and  at  their  neighhoring  points,  hhd^ 
accordingly  is  defined  only  for  such  point's  of  the  region,:  which 

are  also  interior  points  of  the  region  (see  (5))i  the  "tothiiiy 

of  such  points  ve  designate  hy  .  We  then  obtain  dreen's  fomiia 

k*  Au-Av  s  y*R(M),  (6) 


Vk 


where  R  (M,)  is  a  defltiable  ilhsar  difference  expression  for  each: 
point  of  the  houndary.,  strip  ^1^  +  »  thou^i  we  omit  its  more  ; 

precise  form  i  ■  f**/  indicates  in  this  .case  the  set  of  boundary  points 

I  a.  Boundary  Value  problems  and  Eigenvalue  Problems  V 
1*  jSj*  Theory  of  ]fem«dary  Value  Problem,  the  boundary 
value  problem  ibr  llneai*  elliptic  homogSnaous  difference  equations  of 
second  order,'''whieh.'  correa^ads;'  ,j^’'^he'’hias6i(Jal-  bouadRry^  V 
for  partial  diff^imniial;  equations,. 'bah  be'  ;fafmwlatad-ai^.vfj^^ 
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In  a  grid  region  (S  let  there  he  given  a  self-adjoint 
elliptic  linear  diffei^aoe  expression  of  the  second  order  L  (14). 

It  nay  originate  from  a  quadratic  expression  B  (H ,  14)  which  is 
positive  definite  in  the  sense  that  it  cannot  vanish  if  U.  and 
Wy  themselves  do  not  vanish. 

We  now  determine  in  Gj,  a  function  tC  satisfying  the  difference 

** 

equation 

L(lt)  0  , 

and  which  coincides  with  prescribed  values  at  the  boundary  points  of 
this  grid  region. 

Our  requirement  is  representad  by  just  as  many  linear  equations 
and  therefore  function  values  that  are  to  be  determined,  as  there  are 
interior  grid  points  of  the  grid  region.^  Some  of  these  equations  are 
honogeneoua,  namely,  those  which  Correspond,  to  grid  points  which,  with 
their  four  neighbors,  lie  in  the  interior;  others,  in  which  boundary 
points  Of  the  grid  region  enter,  are  nonhomogeneoue .  if  we  set  the 
right-hand  side  of  the  nonhomogeneous  system  of  equations  ( 1 .e . ,  the 
boundary  valuef  of  U)  equsl  to  zero,  then  it  follows  at  once  from 
Green’s  formula  (1)  that  B  (W ,  vanishes  when  w^^  set  and 

because  of  the  positive  definite  character  of  there  also 

follows  the  vanishing  of  liff,  and  therefore  also  of  14.  Thus  the 
difference  equation  has  the  solution  K®®  if  the  boundary  values 

v^ish;  or,  in  ether  words,  the  solution  is  xmlquely  deteralnod,  If  | 

^ ^ 

**  ,  3^  considers  an  arbitrary  difference  e<iuatlon  of  the  8e(f0nd 

order  lU  ( 14 )  s  0  as  a  system  of  linear  equations,  and  conB;bruct,9  thi 
trahepoeed  system  of  equations,  then  this  transposed  set  is  represented 
by  the  adjoint  difference,  equation  M  (V) «  O  .  Thus  the  above  self- 
adjoint  difference  equation  gives  rise  to  a  linear  eipsatlon  system  with 
a  ■symmetric  matr’ix. 
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at  all,  by  the  boundary  values,  since  the  difference  of  two  solutions 
with  the  same  boundary  values  must  vanish.  If,  however,  a  linear 
system  of  equations  with  Just  as  many  imknowns  as  equations  possesses 
the  property  that  for  vanishing  right-hand  sides  the  unknowns  them¬ 
selves  must  also  vanish,  then  the  f^lndamental  theorem  of  the  theory 
of  equations  asserts  that  for  an  arbitrarily  preaci'-ibed  right-hand 
side  exactly  one  solution  must  exist.  In  our  case  there  follows  at 
once  the  existence  of  a  solution  of  the  boundary  value  problem. 

Ws  see  that  for  our  elliptic  difference  equaticus  tne  vinique 
determination  and  the  existence  of  a  solution  of  the  boundary  value 
problem  are  related  to  one  another  by  the  fundamental  theorem  of  the 
theory  of  Ilneai*  equations*,  while  In  the  theory  of  partial  dlfferen- 

V 

tla^  equations  both  facts  must  bo  proved  by  quite  different  methods. 

The  basis  for  this  difficulty  is  to  be  found  in  the  feet  that  the 
differential  equations  are  no  longer  equivalent  to  o  finite  number 
of  equations,  and  so  one  can  no  longer  depend  upon  the  equality^of 
the  number  of  imknovns  and  equations. 

81nce  the  difference  equation 

"  A  t*.  ss  0  ■ 

originates  from  the  positive  definite  quadratic  expression 

k‘  *  «) )  . 

then  the  boundary  value  problem  of  this  dlffe  cneo  equation  is  always 
xiniquely  solvable  . 

The  theory  for  difference  eq..abionB  of  higher  order,  e.g.,  of 
fourth  order,  is  developed  entirely  >  .\rallel  to  that  for  the  difference 


enuations  of"  second  order^  for  which  tho  example  of  the  difference 


equation 


AA'U  --  o 


may  be  sufficient.  In  this  case  the  values  of  the  function  14  in  the 
boundary  strip  must  be  given.  Evidently  here  also  the 

difference  equation  yields  .just  as  many  linear  equations 

as  there  are  unknovn  function  values  at  the  points  of  In  order 

to  demonstrate  the  uniqjieness  of  the  solution  of  the  boundary  value 
problem  ve  need  only  show  that  a  solution  which  has  the  value  aero  in 
the  boundary  atrip  -f  necessarily  vanishes  identically.  For 
this  purpose  we  note  that  the  sum  over  the  corresponding  quadratic 


eiiqpresBlon 


(Alt)' 


for  such  a  function  vanlshesi  as  one  laay  notice  by  transforming  the 
sun  aooordlng  to  Oreen's  fonaula  (6).  Vhe  vanishing  of  the  sum  (7)> 
however^  causal  the  vahlshlng  of  in  all  points  ,  and  according 
to  the  above  proof  this  can  only  happen  for  vanishing  bo^dary  values 
if  the  function  1i  assumes  the  value  zero  throughout  the  region.  Thus« 
our  assertion  is  proven  and  the  uniqueness  of  the  solution  of  the  boundary 
value  problem  of  the  difference  expression  is  guaranteed. 5 


5  For  a  discussion  of  the  solution  of  our  boundary  value  problem 
by  Iterative  methods,  see,  among  others,  the  treatment  in:  "uber  Rand- 
wertaufgaben  bel  partiellen  Dlfferenzenglelchungen",  by  H.  Courant, 
Zeitsehr.  Migaw,  Mathenatlk  u^.  MacliatJ.k,  6  (1925)^  PP-  322-325.  in 
addition  one  is  referred  to  the  report  of  E.^Henky,  Zeitsehr.  f.  angew. 
Math,  u.  Mach..  2  (ISSS),  p.  58  ff. 
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2.  Relatioaa  to  Mlnlaum  Problema .  The  above  boxmilary  value 
problem  is  related  to  the  following  minimum  problem:  among  all 
functions  y)  defined  in  the  grid  region  which  assume 

prescribed  values  at  the  boundary  points,  there  is  to  be  found  such 
a  function  ,  y  )  for  which  the  sum  over  the  grid  region 

assumes  the  smallest  possible  value.  In  this  case  we  assume  that  the 
quadratic  difference  expression  of  the  first  order  is  positive 

definite  in  the  above  sense.  (Sea  p.  6.)  One  may  notice  that  the 
difference  eqjaation  L(<e)s0  results  from  tnls  minimum  req[uirement 
as  a  restriction  on  the  i^olhtton  (jps  “*<  (  X  ,  y  ),  where  L  (  the 

difference  expression  of  the  sebond  order  derived  in  the  above  manner 
from  either  by  applying  the  rules  of  the  differential  cal¬ 
culus  to  the  sum  as  a  function  of  a  finite  number 

of  values  of  ^  at  the  grid  points,  hr  elmllarly,  by  employing  the 
usual  methods  from  the  calculus  of  variations , 

By  way  of  example,  the  boimdary  value  problem  of  finding  a 

II 

solution  of  the  eqpatlon  O  which  assumes  prescribed  boundary 

values  is  eqpivalsnt  to  the  problem  of  minimizing  the  sum 
over  all  those  functions  which  assume  these  bomdary  values. 

The  situation  is  entirely  similar  for  difference  equations  of  the  ' 
fourth  order,  where  again  we  restrict  ourselves  to  the  example  of 
AAV’aO  •  The  boundary  value  problem  corresponding  to  this  difference 
equation  Is  equivalent  to  the  problem  of  minlmlBlng  the  sum  K’SS  dv) 

h 
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over  all.  thoee  functions  W(  X,  which  take  on  the  presci:!  cc'  vaLur-B 
in  the  boundary  strip  .  Besides  this  sius,  still  other  ex-prei-i.-.ilons 
which  are  quadratic  in  the  second  derivatives  yield  the  equation 
6  Al<.  =  o  under  the  requirement  that  they  are  to  be  minimized,  as 
for  example  the  sum 


in  which  there  appear  second  difference  quotients  exclusively  at 
points  of  .  . 

That  the  poe^^d  minlinum  problem  always  possesses  a  solution 
follows  from  the  theorem  that  a  continuous  function  of  a  finite  number 
of  variables  (the  functional  values  of  at  the  grid  points)  must 
always  have  a  mlntuum  If  this  function  is  bounded  from  below,  and  If 
It  tends  to  infinity  ««  soon  as  at  least  one  of  the  independent 
variables  tends  to  infinity,^ 

3.  The  green's  Rmotion.  One  can  treat  tha  boundary  value, 
problem  of  the  nonhomo^neous  equation  L  ( U ~  f  similsfly  to  the 
boundary  value  problem  of  the  homoi^naouB  equation  L.  (I(  )s  0  .  Xt  is 
sufficient  for  the  case  of  the  nonhomogeneous  equation  to  liaii^  oneself 
to  the  ease  when  the  boundary  ^nilues  of  ti  vanish  everywhere,  since  we 
obtain  the  solution  for  other  boundary  values  by  the  addition  of  a 
suitable  solution  of  the  homogeneous  equation.  To  solve  the  system  of 
linear  equations  which  is  represented  by  the  boundary  value  problem  of 
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|_  (U)k  -f  ,  we  first  choose  the  function  -f  (  X  >  y)  so  that  at  a  grid 
point  with  the  coordinates  X  s:  ^  ^  ~  ^  It  assumes  the  value  -  n  , 

and  at  all  other  grid  points  assumes  the  value  zero.  If  K (  ^ 5  ^  >  *1 ) 

is  a  solution  of  the  special  resulting  difference  equation  which  will 
now  depend  on  the  parametric  point  ( ^  H  which  vanishes  at  the 
boundary,  then  the  solution  corresponding  to  an  arbitrary  function  is 
represented  by  the  sum 

Wo  caU  the  function  K  (  Xu  y  i  |  #  )>  depends  upon  the 

points  (X  ,y)  and  ($  ,>l),  the  arson ’s  function  of  the  difference 
expression  L  ( u) .  If  we  denote  by  fc  (  X  ,  y  i  ?  » »l )  the  Green » s  function 
of  the  adjoint  eaipreaslon  M  ('»’ )>  'tben  we  have  the  relation 

Which  can  be  obtained  iaoediately  from  Green's  formula  (5)  If  one  sets 
*1)  VsKi  X,yj5  ,►[).  J^r  a  self-adjoint 
difference  expression  there  results  froa  the  above  relation  the  symmetry 


relation 


K  (?,![;  5.1) 


K(s,n)  in ) . 


4,  tlaeavalue  I>robleas.  Self -adjoint  difference  e3cpr«»*ion8 
Leu)  give  rise  to  eigenvalue  problems  of  the  following  types  let 
there  be  sought  the  values  of  a  paraseeter  \  --  the  eigenvalues  —  for 
which  the  difference  equation  in  (»j^ 
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possesses  a  solution  --  the  eiganfunctlon  —  vhlch  vanishes  on  the 
boundary  H  ^  . 

The  eigenvalue  problem  is  eq,uivalent  to  the  principal  axis 
problem  of  the  quadratic  form  BCMi'U)*  There  are  just  as  many  eigen¬ 
values  as  interior  grid,  points  in  the  region  , 

and  Just  as  many  corresponding  eigenfunctions  **'  i  The 

system  of  eigenfunctions  and  eigenvalues  emd  their  existence  results 


from  the  minimum  problem; 

Among  all  the  functions  (^( y>  ^ )  which  vanish  on  the  boundary 
and  satisfy  the  (in- i  )  orthogonality  conditions 

and  the  normalization  condition 


h‘£-£/  = 


we  seek;  that  one  <(>  ■  li  for  which  the  sum 

k*  sz 

<3(, 

assumes  the  smallest  value .  The  value  of  the  minlinuo  is  the  eigen¬ 
value  and  the  function  for  which  It  is  asaianed  is  the  eigenfunction. 

Because  of  the  orthogonality  of  the  eigenfunctions,  i.e., 
k  every  grid  function  9  (  X  ,  y  .)  wh,ich  vanishes 

on  thir  boundary  can  be  developed  in  a  series  of  eigenfunctions  In  the  form 

where  the  coefficients  ^  obtained  from  the  equation 

®h 

m  this  way  we  obtain  in  particular  the  following  representation 
of  the  Green's  function: 

It  JktO) 
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§  5.  Relation  to  the  Problem  of  a  Random  Walk® 

Our  topic  is  related  to  the  question  of  probability  calcu¬ 
lation,  namely  the  problem  of  a  random  walk  in  a  bounded  region.^ 

In  a  grid  region  one  may  Imagine  the  grid  lines  as  paths  along 

which  a  particle  can  move  from  a  point  to  a  neighboring  point.  In 
this  path  -  network  our  particle  may  now  move  at  random,  choosing  by 
chance  at  each  path  junction  among  the  four  available  directions  — 
all  four  are  equally  probable.  The  random  walk  stops  as  soon  as  a 
boundary  point  of  la  reached,  where  the  particle  may  be  absorbed. 

We  shall  ask: 

a  ‘ 

1.  What  Is  the  probability  V^( P >  R )  that  starting  from  a 
point  P  one  can  arrive  in  some  time  or  other  by  a  random  path  at  a 
boundary  point  R  ? 

2.  What  is  the  mathematical  expectation  V (  PjQ)  that  by 
such  a  random  path  starting  from  P  one  reaches  a  point  Q  of 
without  first  meeting  the  boundary? 

We  may  consider  this  probability  or  the  mathematical  expectation 
more  precisely  by  the  following  process.  We  imagine  a  unit  quantity 
of  some  Substance  existing  at  a  point  P  .  The  substance  may  spread 
out  in  our  path  net  with  a  constant  velocity,  travelling  by  chance  one 
grid  space  in  one  time  unit.  At  each  grid  point  exactly  a  quarter  of 

® Here  {  5  is  imneaeasary  for  the  discussion  of  the  limiting 
process  In  f  4. 

9  There  is  an  essential  difference  in  the  way  in  which  the  boun¬ 
daries  of  the  region  are  introduced  in  the  following  consideration  and 
in  well-lmown  methods,  such  as,  for  example,  those  which  have  bean  used 
in  connection  with  Brownian  molecular  motion. 


i 


I 
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ths  incoming  substance  sbculfi.  flow  toward  each  of  the  four  directions. 
The  amount  of  substance  which  arrives  at  a  boundary  point  will  be  held 
there.  If  the  source  point  P  is  a  boundary  point,  then  the  entire 
quantity  of  substance  will  remain  there. 

We  Interpret  in  general  the  probability  w(  P ;  fl)  of  reaching 
the  buuHdni-y  point  ft  by  a  landom  walk  starting  from  F  without  having 
before  touched  the  boundary,  as  the  quantity  of  substance  which  collects 
after  an  Infinite  tijiie  at  this  boundary  point. 

The  probability  (  P;^?)  of  reaching  a  point  Q  In  exactly 
n  steps  from  point  P  without  touching  the  boundary,  we  interpret 
as  the  quantity  of  substance  at  point  Cj  after  tt  time  units,  when  P 
and  are  both  Interior  points;  if  P  or  is  a  boundary  point,  then 
ve  set  it  equal  to  zero. 

The  value  of„Ej,(P;<?)  ia  just  the  amount  moving  in  M  stepa  , 
from  P  toward  Q  without  meeting  the  boundary,  divided  by  4;  there¬ 
fore  it  follows  that  E  (  P  ;Q  )  ar  6^(C5  j  P  ). 

By  the  mathematical  ex^ctatlon  v(P;(?)  for  the  above  mentioned 
random  walk  going  once  from  P  to  point  ^  ,  we  mean  in  general  the 
infinite  sura  of  all  these  probabilities,^® 

p*  0 

which  is  therefore  for  interior  points  P  and  (p  the  sum  of  all  quanti¬ 
ties  of  substance  which  have  pasaed  the  point  CJ  at  different  moments 


ro 


We  shall  prove  their  convergence  shortly. 
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of  time.  Therefore  the  expectation  value  1  will  he  assi.gnecl  for  reach 
ing  point  .  For  boundary  values  this  expectation  is  equal  to  7.ero, 
Denoting  the  quantity  arriving  at  a  boundary  point  R  in 
exactly  n  steps  by  Fj^(  P;  R),  the  probability  W'(  P,  R)  is  thus 
represented,  by  the  infinite  sum 

w(e;R)  =  ^  Fi)(f;K)  , 

i;>wo 

all  of  whose  members  are  positive,  and  whose  partial  sums  can  never 
be  greater  than  one,  because  the  substance  arriving  at  the  boundary 
will  make  up  only  a  part  of  the  original  quantity  of  substance.  But 
with  this  the  convergence  of  the  series  is  guaranteed.  . 

One  can  easily  see  that  the  probability  P  ;  Q),  i.e.,  the 
quantity  of  substance  arriving  at  a  point  <5  after  exactly  W  steps, 
tends  towai’d  zero  with  increasing  H  ,  itor  if  Q  is  a  point  from  which 
a  boundary  point  R  will  be  reached  in  exactly  Hi  steps,  and  if  we 

have  E^(PjQ)>e(>0  ,  then  after  in  steps  at  least  the  amount 

**IVI 

arrives  at  this  boundary  point  R  j  but,  because  of  the  conver¬ 
se 

gence  of  the  sum  (f  ;ll)»  the  quantity  of  substance 

arriving  at  the  boundary  point  n  tends  toward  zero  with  time,  and 
the  value  of  E^(  P;<J)  Itself  must  therefore  also  tend  toward  zero 
with  Increasing  n  ;  i.e.,  the  probability  of  an  Infinitely  long  path 
remaining  in  the  interior  is  zero. 

From  this  it  foi:i.ovs  that  the  whole  quantity  of  substance  must 
finally  arrive  at  the  boundary)  in  other  words,  the  sum  extended  over 
all  boundary  points  R  is 

S  v-C  P;  R)  =  I  . 
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We  still  hRVR  to  prove  the  convergence  of  the  Infinite  series 
■for  mathematical  expectation  v(  PjQ), 

v(P;<?)=^  Si  E^rP;<?), 

i)si> 

For  this  purpose  wa  notice  that  the  quantities  £^(  P  ) 
sattsf^r  the  following  relation 

+  6HCf;9,)  +  6,(P;9^)|  , 

where  to  are  the  four  neighboring  points  of  the  interior 
point  <9  .  That  is,  the  quantity  of  substance  arriving  after  lt+( 
steps  at  point  (?  is  one  fourth  of  the  amount  of  substance  reaching 
the  four  neighboring  points  of  Q  after  W  steps.  If  one  of  the 
neighboring  polnta  of  Q  is  a  boundary  point,  e.g.,  R  #  then 

there  follows  that  no  quantity  of  sutistanoe  flow#  from  this  boundary 
point  to  <J,  since  in  the  expression  we  have  set  E^CPjR)  equal  to 
zero.  Purtheimore  for  an  interior  point  E^(  P;  P)*  I  ^  and.  of 
course 

From  these  relations  we  obtain  for  the  partial  sums 

OssO 

the  equations 
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if  P  does  not  coincide  with  Q  ;  otherwise  we  have 

V,*,(P;P)=  I+^K(P;P,)  +  v„(P;Pj)+v;CPjP,)  +  Vh(P;P^)|  , 

i.e.,  the  expectation  of  a  particle  coining  back  to  its  starting  point 
is  composed  of  the  expectation  of  reaching  the  point  P  again  on  a  non¬ 
disappearing  path,  namely  CP-J  ?4)|, 

and  from  the  expectation:  unity,  which  eaqipesBes  the  fact  that  the  whole 
substance  was  originally  present  at  this  point. 

Therefore  the  quantities  VJ^CPiQ)  satisfy  the  following 
difference  equations  ir'^ 

A  ( P>  (?  )  »  ±  E  „  (  P;  <} )  ,  it  P  #  (J  , 

h  " ... 

■  ■  ^  \  ..  , 

Here  V*^  (  p  j  <} )  is  eqjial  to  sero  if  Q  is  a  boundary  point. 


In  this 'base  the  ^-operation  is  related  to  the  variable  point  Q  . 

This  eqjiition  may  be  interpreted  as  a  type  of  heat  conduction 
equation.  Namely,  if  one  considers  the  function  Va(PiQ)  »  function 
of  time  "t  which  is  proportional  to  n  so  that  "t  a  n't  and  Vk(P»Q) 
a  V'(P;<9;'t  )a  vCt  ),  instead  of  as  a  function  of  the  index  n  as  in 
the  above  representation,  we  can  then  vrlte  the  above  equations  in  the 
following  form: 

a,v(t)  =  is.  p+Q, 

h* 

avC*)  »  ^  ^  p=tj 

For  the  limiting  process  of  a  similar  difference  equation  to  a  parabolic 
differential  equation  see  Part  II,  §6,  p.  55* 
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The  solutiOii  of  this  boundary  value  problem  is,  as  already 
demonstrated  earlier,  uniquely  defined  for  any  right-hand  side  (see 
p.  10)  and  depends  continuously  on  the  right-hand  side.  Nov  since 
the  values  of  E  P j  9)  tend  toward  zero,  therefore  the  solutions 

V^CPjQ)  converge  toward  the  solutions '\r(  P ;  (J )  of  the  difference 
equations 

4v(P;(?)  =  0  ,  it  , 

avCP;0)  =  "^1  '  P  “  ^  > 

with  the  boundary  values  V(P;R)*0. 

Therefore  we  see  that  the  mathematioal  expectation  lr( P  j  Q) 

.  .« 

exists  and  is  nothing  more  than  the  Green's  function  K(Pi9)  "belong¬ 
ing  to  the  difference  equation  but  supplied  with  the  factor 

4.  The  symmetry  of  the  Green's  function  K(  P)9)~  IC(9^  is  an 
immediate  consequence  of  the  symmetry  of  the  values  pjQ  ),  with 
whose  help  it  was  defined. 

The  probability 'Vtf(  P  J  R)  satisfies  the  relation 

'W'(P;  R)«  ^|tr(P/;R)+w(P»;RHw(P,;a)4V‘Cp4;R)| 

with  respect  to  p,  and  therefore  also  the  difference  equation 

A'Ur  SB  0 

Since  Pj  ,  >  ^3  >  ^4  neighboring  points  of  the 

interior  point  P,  therefore  each  path  from  P  to  R  must  pass  through 
one  of  these  four  point;?,  and  each  of  the  four  path  directions  is 
equally  probable.  Further,  the  probability  of  going  from  one  boimdary 
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point  K  to  fu'iothsr  R  Is  zero,  “Ur  ( R  ,  fZ  )*=0  ,  and  if  the  two  points 
R  and  R^cOincide,  we  have  w(R,R)«  I  .  Therefore  'UrCPjR)  is  the 
solution  Of  the  hoiuidary  value  problem  Avr»0,  where  the  boTindary 
value  1  is  prescribed  at  the  boundary  point  R  and.  the  value  0  is 
prescribed  at  all  other  points  of  the  boundary.  The  solution  of  the 
boundary  value  problem  for  arbitrary  boundary  values  I4.  (  R  )  then  has 
simply  the  form  'U,(  P  )«^vr(  P  ;  R)  K  (  R),  which  is  to  be  stunmed 
over  all  boundary  points  R  If  we  here  substitute  for  It  the  function 
'Us  I  ,  we  thus  again  obtain  the  relation  \"  S  W’(  PjR). 

The  present  Interpretation  of  Green's  function  as  expectation, 
permits  further  properties  to  be  Immediately  recognized.  We  mention 
only  the  fact  that  the  Green's  function  increases  if  one  changes  from 
the  region  G  to  one  which  contains  G  as  a  partial  region;  for  then 
the  number  fv  of  the  possible  grid  paths  leading  from  one  point  P  to 
another  Q  without  touching  the  boundary  also  increases. 

Naturally  for  more  than  two  independent  variables  corresponding 
relations  hold.  We  satisfy  ourselves  with  the  observation  that  other 
elliptic  difference  equations  permit  a  similar  probability  interprota" 
tion. 

If  one  goes  through  the  limiting  process  to  vanishing  meph  width, 
which  is  easy  to  perform  by  the  methods  of  the  following  paragraph, 
then  Green's  function  for  the  grid  changes  to  within  a  nufflerleal  factor 

Tp'  '  ~  •  '  "  ■■■  ■--■  ^  -n  — ■  ■  1  -  —  -  Ml- 

One  can  easily  see  that  the  probability 'Kr( P ; R )  of  reaching 
the  boundary  is  the.  boundary  expression  R  (K(  PjQ))  written  with  re¬ 
spect  to  Q  with  the  Green's  function  K  (  P ;  Q ),  by  identifying  ■»<.(  X  ,  ij  ) 
in  Green's  formula  (5)  with  v(  Pj  Q)  and  v(  X  ^  y  )  with  v(  P ,  Q). 
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to  the  Green's  function  of  the  potential  equation;  a  similar  relation 
exists  between  the  expression  fi  *'W{**;  R)  and  the  normal  derivative  of 
Green's  function  at  the  boundary  of  the  region.  In  this  way  the 
Green's  function  of  the  potential  equation,  for  example,  could  be 
Interpreted  as  the  specific  mathematical  expectation  of  going  from 
one  point  to  another  without  touching  the  boundary. ^5 

After  the  limiting  process  from  the  grid  to  the  continuum,  the 
influence  of  the  prescribed  grid  directions  of  the  random  walk  vanishes . 
This  fact  is  of  Ij^portance  when  one  undertakes  the  limiting  process 
with  a  more  general  random  walk  problem  without  prescribed  directions, 
and  is  in  principle  sm  Interesting  problem;  however  it  exceeds  the 
scope  of  this  disousslon,  but  is  one  to  which  we  hope  to  return  on 
another  occasion. 

§  Limiting  Process  to  the  Solution  of  the  Differential  Eqjiation 
1.  The  Boundary  Value  Problem  of  Potential  Theory.  In  carrying 
out  the  limiting  process  for  the  solution  of  the  difference  equation 
problems  to  the  solution  of  the  corresponding  differential  equations, 
we  will  not  seek  the  greatest  possible  generality  in  the  formulation 
with  respect  to  the  boundary  and  the  boundary  values  themselves ,  in 
order  to  make  the  characteristics  of  our  method  clearer. Accordingly 


I? 


Thereby  the  area  of  a  surface  is  assigned  as  the  expectation  of 


reaching  the  surface  element. 
l4 

It  might  be  remarked  that  the  extension  of  our  method  to  more 
general  boundaries  and  prescribed  boundary  values  in  no  way  gives  rise 
to  fundamental  difficulties. 


I 


we  assume  that  a  simply  connected  region  G  is  given  in  the  plane, 
whose  boundary  consists  of  a  finite  number  of  ares  having  contlmio\xs 
tangents.  In  a  region  containing  G  i^^  its  interior  let  there  be 
given  a  continuous  function  -f  (x  >  y  )  having  continuous  partial 
derivatives  of  fl;*Gt  and  second  order.  Jbr  the  grid  domain  , 

•  #  j"ii 

corresponding  to  the  mesh  width  fl  and  to  the  region  to ,  let  the 
boundary  value  problem  of  the  difference  equation  be 

solved  with  the  same  boundary  values  which  are  assumed  by  the 
function  -f  (  X  >  y )  at  the  boundary  points  of  (a^>  call  this  solu¬ 
tion  ttj^(X»y).  We  wish  to  prove  that  with  vanishing  mesh  width 
K  the  grid  function  IC  converges  toward  the  solution  1i  of  the 
boundary  value  problem  of  the  partial  differential  equation 

for  the  region  Q ,  where  the  boundary  values 
for  the  region  G  arc  again  provided  by  those  values  which  the 
function  •?  (  X ,  y)  assiames  on  the  boundary  of  G  .  Furthermore  we 
shall  show  that  for  every  region  lying  entirely  In  the  Interior  of 
G  the  difference  quotients  of  14. of  arbitrary  order  converge 
uniformly  toward  the  corresponding  partial  differential  quotients 
of  the  limiting  function  14  (  K,  y  ). 

In  carrying  through  the  proof  of  convergence  it  is  convenient 
to  replace  the  requirement  that  ‘U(x  >  y  )  assumes  the  boundary  values 
by  the  following  weaker  condition:  If  Sj.  is  that  boundary  strip  of 
th*'  region  G  whose  points  are  at  a  distance  Iccc  than  {"  from  the 
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botmdary,  then  the  integral 

I 


sV 

tends  toward  zero  with  decreasing  f 

Our  proof  of  convergence  depends  upon  the  fact  that  for  every 
subregion  Q*  lying  entirely  within  the  Interior  of  the  region  G  , 
tha  function  • » *J  )  *hid  every  difference  quotient  remains  bounded 

with  decreasing  h  and  are  "eqtii-continuous"  in  the  following  sense: 
There  exists  for  each  of  these  functions  V  j^(  X,  ^  )  a  qpemtlty  3  (6  ) 
dependent  only  on  the  region  and.  not  on  such  that 

I  ^  > 

when  the  grid  pp,lnt  P -and  P,  of  the  ^id  region  6^  He  In  the  given 

•  '  il  '  *  .  ^ 

subregion,  and  are  at  a  distance  less  than  S  {€)  from  one  another. 


Once  we  have  proved  the  asserted  equi-continulty  then  we  can 
select  by  known  methods  a  subse’Tuenoe  of  our  functions  in  such 
8  way  that,  together  with  its  difference  quotients  of  every  order,  it 
converges  uniformly  in  every  subregion  Q*  toward  the  Halting  function 


^  Note  that  our  weaker  bovaidary  Value  requirement  actually  provides 
the  unique  characterisation  of  the  solution,  as  follows  from  a  theorem 
which  is  easy  to  prove;  If,  for  a  function  satisfying  the  differential 
equation  O  in  the  interior  of  S  ,  the  above  form  of 

the  boundary  condition  is  satisfied  with  X,  ^)«0,  and  if 

5f((  ?^/9x)*4  exists,  then  ti( X , ^  )  is  identically  zero. 

(See  Courant,  "ih»r  die  liSsungen  der  Dlff ,  01.  der  Physlk",  Math,  Annalen. 
85.  especially  pp.  296  ff.) 

In  the  case  of  two  independent  variables  the  fact  that  these  bound¬ 
ary  values  are  actually  taken  on  can  be  concluded  from  our  weaker  re<^ire- 
mentj  in  the  case  of  more  variables  one  cannot  expect  the  corresponding 
result  in  general,  because  on  the  boundary  there  may  be  exceptional  points 
at. which  boundary  values  need  not  be  assumed,  while  for  the  weaker  require 
ment  a  solution  still  exists. 
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or  its  differential  quotients,  reBpectlvely.  The  limiting 
function  possesses  corresponding  derivatives  of  arbitrarily  high  order 
iT'  each  interior  subregion  Q*  of  $  ,  and  there  satisfies  the  partial 
differential  equation  ?*tt  /^X*  +  s=0  .  If  ve  then  show  that 

it  Batlsfies  the  boundary  conditions,  we  recognise  it  in  the  solution 
of  o\ar  boundary  value  problem  for  the  region  .  Since  this  solution 
is  uniqtiely  determined  then  it  is  farther  shown  that  not  only  a  sub¬ 
sequence  of  the  functions  but  that  this  sequence  of  functions 
itself  also  possesses  the  asserted  convergenoe  property, 

The  equl-eontinuity  of  our  quantities  follows  from  the  demon¬ 


stration  of  the  following  results: 

1,  As  h  decreases  the  sums 

extended  over  the  grid  regiOD.  are  bounded. 

S,  If  satisfies  the  difference  equation  AwaiO  In  _ 

a  grid  region  and  if  with  decreasing  h  the  sum 

h*2Z  w*  , 

_  ■ 

extended  over  a  grid  region  (5^  eorrespondlng  to  a  subregion  Q* 
of  G  is  bounded,  then  for  every  fixed  subregion  lying  entirely 
in  the  interior  of  (J*  the  sum 


Here  and  when  later  convenient  we  shall  drop  the  subscript  h 
from  grid  functions. 


a^uj  liklit  iiiillk  ill.d^ 


I- 


extended  ovi.r  t,-c  corrcEnondt-.ij  region  remains  bounded 

K 

vith  decreasing  K  . 

Together  vith  1,  it  follows  from  this  that,  since  all  of  the 
difference  quotients  W  of  the  function  M  ^  satisfy  the  difference 
equation  AwaO,  each  of  the  sums 

I  k*  2Z  -w* 

s* 

is  bounded,  n 

3*  Prom  the  boundedness  of  these  sums  there  follows  finally 
the  boundedness  and  equi-continuity  of  all  the  difference  quotients 
themselves. 


2.  Proof  of  the  LemmaB.  The  proof  of  theorem  1,  follows  from 
the  fact  that  the  function  values  are  themselves  bounded.  Since 
the  largest  and  smallest  value  of  the  function  will  be  assumed  on  the 
boundary^T,  it  therefore  tends  toward  preecribed  finite  values,  The 
boundedness  of  the  sum  is  an  ismediate  result 

of  the  alaimuD  property  of  our  grid  Auiction  formulated  in  $  2,  2, 
whereby 

certainly  holds.  The  sum  on  the  right  converges  vith  decreasing  mesh 


'  Wo  remarlt  ejsplicitly  that  vith  regard  to  the  application  of  the 
method  to  other  differential  aq^atlon8,  we  need  not  use  this  condition. 

In  this  connection  we  could  use  the  inequality  (15)  or  employ  an  alterjia- 
tlve  method.  (Translator’s  note;  see  Courant  and  Hilbert,  Methoden  der 
mathematlschen  Phyalh.  Vol.  1,  Chapt.  ill,  J,  "Alternative  to  the  Theory 
of  Integral  Equations"). 
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width  toward  the  Integral 
aasuraed  to  exist. 


which  we  have 


In  order  to  prove  the  lemmas  formulated  under  2,  we  consider 
the  quadratic  sum 

Mfi 

where  the  sum  is  to  extend  over  all  interior  points  of  a  square  j 
(see  Fig.  1).  We  designate  the  function  values  on  the  outer  boundary 
of  the  square  Q,  by  fer,',  those  on  the  second  boundary  line  S,  by  W^. 
Then 


Fig.  1. 

Green ' s  ,  formula  yields 

h’S!  (wj+wjtw'+w,  )=  £Cw’-w-.*)S£w’ -Sw?,  (8) 

0,  s  s,  s. 

where  the  sunnaation  on  the  right  is  to  extend  over  the  two  outer 
boundary  lines  and  where  "V,  and  'Wjj  refer  to  neighboring 

points.  We  now  consider  a  set  of  concentric  squares  '  Qx  * 

•  • •  >  Q u  with  boundaries  >  S  ,  $  each  of  which  arises 


1 

I 
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from  the  previous  one  in  such  a  way  that  the  bordei'  of  nearest 
neighboring. points  Is  added.  (See  Fig.  1.)  To  each  of  these  squares 
we  apply  the  appraisal  (8)  and  notice  that 

k*2SCw,%w5+>»,‘+W5’') 

’  % 


X  ^  "'s  ^ 

holds  for  jC  ^  !  »  If  ve  add  the  series  of  h  tnegualitles 


<?• 


(o^kcn)  , 


we  obtain 


ZE  (wicVH^)  ^  Z  w*  -  Z  w-*  «  s  V* 

Q.  s,, 

This  Inegyiallty  we  stun  from  I  to  ns  N  •  One  thus  obtains 


N*h' 


<»• 


w 


In  which  the  sum  bh  the  right  extends  over  the  entire  sqjuare  , 

Qy  reduction  tf  the  mesh  width  we  now  let  the  squares  and 
Ql^  tend  toward  two  fixed  aonoentrie  squai^ea  lying  In  the  interior 
df  (9  and  a  dlatanoe  <t  apart,  ao  that  NK  converges  toward  a  and 
We  find  that  " 

"V  ■:  (9) 

■  I?.  "  '*  Qm 

holds  Independently  of  the  mesh  width, 

This  inequality  —  for  sufficiently  small  mesh  width  —  holds 
of  course  not  only  for  the  two  squares  and  but,  with 

another  constant  CL  ,  for  any  two  subregions  of  ig  such  that  one  is 


W 

..uV. 
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contevined  entirely  within  the  other.  Therefore  the  theorem  of  2.  is 

1  Q  ' 

proved. 

In  order  to  prove  the  third  result,  that  In  each  interior  sub- 


region  the  function  tx.^»nd  all  of  its  difference  qjiotients  ^ 
remain  bounded  and  equi-continuous  with  the  refinement  of  the  mesh 
width,  we  consider  a  rectangle  R  with  comer  points  p®  ,  Qs  ,  P  ,  Q 
(see  whose  sides  Oo  ®hid  PO  sre  psrnllel  to  the  X  -axis 


and  have  the  length  CL  . 


Fig.  2 

We  start  with  the  relation 

w{(?.)--w(p.) » h  2  +  K*  ss;  ■w,, 

f*  R 

and  with  the  inequality  which  follows  directly 

l-wc9o)-v»(p.)l<K2.1’'xl  +  k^2£lw»sl.  (11) 
'  '  n  R 


lb  If  we  do  not  assume  that  A  WsO  >  +'beh  in  place  of  the  lne(juallty 
{ 9 )  we  find 

G**  <5*  G* 

for  suitable  constants  Cj  and  Independent  of  U,  and  where 
lies  entirely  in  the  interior  of  the  region  Q*,  which  in  turn  is  con¬ 
tained  In  the  interior  of  g . 
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We  then  let  the  Bide  PO  of  the  rectangle  vary  between  an  Initial  line 
P,  Q,  a  distance  b  from  P^Qo  and  a  final  line  Pj  a  distance 
2j>  from  RjCJe  ,  and  we  sum  the  corresponding  4.  |  )  inequalities  (ll). 

k 

Thus  we  obtain  the  result 

|w(p.)-w(<?,)|  ll*  S£1‘^kI  +  , 

in  which  we  extend  the  summation  over  the  entire  rectangle  VRjj  ** 

.  From  Schwarz's  inequality  there  follows 

|aArCft,)~vr<:Co)|^  -r  V  vr^’  +Yilb  .  (is) 

Since  by  aasujiiplion  the  aums  which  occur  here  multiplied  by  k* 


remain  hounded,  it  follows  that  the  difference  1 1gf(  )- V^T  (  )  | 

tends  to  zero  with  the  distajice  Cu,  Independently  of  the  mesh  size 
since  we  can  fix  b  for  each  subregion  G*  of  Q  ,  The  eqiii-continulty 
of  W»  in  the  y -direction  is  thus  proved,  Correspondingly  one 
may  prove  the  sane  for  the  H -direction,  and  hence  for  each  interior 
subregion  Q*  of  G  •  The  bouadadness  of  the  function  G* 

follows  finally  from  its  eqiil-eontlnulty  and  from  the  betindedness  of 


K’22w?. 

G*  ^ 

With  this  proof  one  then  assures  the  existence  of  a  subsecjuenoe 

of  functions  *14  which  converges  toward  a  limiting  function  ‘fi(x  ,y  ), 

«  • 

and  which  does  so  imlforoily  together  with  all  its  difference  quotients 
in  the  sense  discussed  above  for  every  interior  subregion  of  G •  I'hlo 
IJjuiting  function  KC  )  has  continuous  partial  differential  quotients 


I 


of  arbitrary  order  throughout  (5  ,  and  there  satitsfies  the  partial 
differential  equation  of  the  potential 

.  9’ it  _  Q 

5,  The  Boundary  Condition .  To  prove  that  the  solution  satis¬ 
fies  the  previously  stated  boundary  condition,  we  next  show  that  for 
every  grid  function  V  we  have  the  inequality 

K*s:£v*j5  Ai-*lt’Z£(’^^v/)  +  Bhk2;v\  (15) 

*f.h  Sf,h  n. 

where  dcsignntca  that  port  of  the  grid  region  Q  which  lies 

inside  a  boundary  strip  S^'  boundary  strip  (®®®  P*  25) 

consists  of  a.11  points  of  G  whose  distance  from  the  boundary  is  Less 
than  h  j  it  is  confined  between  P  and  a  cui*ve  P|,  Furthermore  A  = 
and  6  are  constants  dependent  only  on  the  region  <5,  and  not  on. the 
mesh  size  h  aor  on  the  function  V . 


r 


Fig.  3.  Fig.  4. 
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In  order  to  prove  the  above  inequality  we  divide  the  boundary 
P  of  G  into  a  finite  munber  of  sections,  for  which  the  angle  between 
the  tangent  and  either  the  X  -axis  or  the  y  -axis  stays  above  some 
positive  value  (such  as  50° )•  In  illustration  let  Jf  be  such  a  section 
of , P  which  is  rather  steeply  inclined  to  the  X-axis,  (See  Fig.  U, ) 
Lines  parallel  to  the  X  -axis  intersecting  the  end  points  of  the 
section  JT  will  cut  a  section  fj,  from  the  neighboring  curve  , 
and  will  define  together  with  J"  and  a  region  of  the  boundary 

strip  That  part  of  the  grid  region  (5  ^  which  is  contained  in 

the  strip  is  called  Sj,  ,  and  the  corresponding  part  of  the 
boundary  P^  is  called  • 

We  now  imagine  that  a  line  parallel  to  the  X  -axis  is  drawn 
through  a  grid  point  of  S^.^.  This  line  will  hit  the  boundary 
at  a  point  P|^  ,  That  particular  section  of  this  line  which  lies 
in  we  shall  call  l*nK  •  Certainly  the  section  will  have  a 

length  smaller  than  CY*  ,  since  T,  is  the  largest  perpendicular  dis¬ 
tance  of  a  point  on  Sf  from  P  .  Thus  the  constant  G  depends  only 
on  the  smallest  angle  which  the  tangent  to  {(■  makes  with  the  X  -axis. 

The  following  relation  holds  between  the  values  of  V  at  the 

i>'(Ph)  =  '^CP|j)  ±  K  SI , 

from,  which,  by  squaring  and  using  Schwarz's  inequality,  one  finds 


point 


and  its  value  at  the  point 


■v-cPi,)  €Zv(F^)  +  2ct'hS!. ■ 
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By  amming  with  respect  to  In  the  )(  -direction,  we  obtain 

KZ  vC\)*+  Zc^r^\i  • 

h 

Suming  once  more  in  the  ^  -direction  yields  the  relation 

5,.  Pt,  Sv 

and  in  order  to  find  the  desired  inequality  (15)  from  this  one  needs 

only  to  set  up  the  corresponding  expression  for  the  other  parts  If  of 
n  ,  and  then  to  add  the  two  express ions. ^9 


We  next  set 


80  that  Vr  .vanishes  on  the  boundary  ,  Then,  since  K  £E(v*f-v*) 

L 

remains  bounded  for  deoreasing  n  ,  we  find  from  (13) 

lv%‘' z  s  s  Kr  ,  (16) 

where  X  l9  a  constant  independent  of  both  the  function  'V  and  the 
mesh  width.  If  ve  extend  the  summation  on  the  left  hand  side  not 
over  the  entire  boundary  strip  'but  rather  over  the  difference 

of  two  such  strips  ,  then  the  inequality  (16)  is  still 

valid  with  the  same  constant  K ,  and  we  can  carry  out  'ths  limiting 
process  for  vanishing  mesh  width.  .'FVoa  the  inequality  (16),  there 


^9  The  same  sort  of  considerations  which  led  to  the  proof  of  the 
inequality  (13)  also  yields  the  inequality 

®K  Ph  ®K 

where  the  constants  C,  and  depend  only  on  the  sregion  Q  and 
not  on  the  mesh  width.  , 


then  results 

■Jr  [f  $  Kr  ,  V-  =  u-f  . 

S>i-5^ 

If  we  now  let  the  smaller  boundary  atrip  approach  the  boundary, 

then  we  find  the  following  Inequality: 

I  fr  fi.  I  «  I  /y  .  I  I  j 

7jj  =  “frjj  cu-t-;  Axdij  ^  Kh  , 

Sr 

which  just  states  that  the  limiting  function  it  satisfies  the  boundary 
conditions  which  we  have  specified. 

4 .  Applicability  of  the  Method  to  other  Problems .  Our  method 
is  essentially  based  on  the  Inequality  (lO)  stated  in  the  above  lemma, 
because  from  it  the  last  two  theorems  of  pages  26-7  follow.  Inequality 
(10)  in  no  way  makes  use  of  special  solutions  or  other  special  proper¬ 
ties  of  our  difference  ejqpresslons,  and  may  therefore  be  Immediately 
03rt;ended  to  the  case  of  arbitrarily  many  Independent  variables,  or  to 
the  eigenvalue  problem  of  the  differential  aquation 
1-  ,  and  we  thereby  obtain  exactly  the  sane  results  concerning 

convergence  properties  as  above. With  some  modifications  which  are 
easily  made  the  method  can  be  extended  to  linear  differential  equations 
of  other  aorta,  particularly  to  those  with  non-constant  coefficients. 


The  esaentlal  distinction  consists  only  in  the  proofs  of  the  boundedness 


Concerning  the  application  of  corresponding  Integral  inequali¬ 
ties,  see  K.  Priedriohs,  "Die  Rand-  und  Eigenwertprobleme  aus  der 
Theorie  der  elastischen  Flatten",  Math.  Annalen,  §8,  p.  222. 


It  is  thus  proved  at  the  same  time  that  each  solution  of  such 
a  differential  equation  has  derivatives  of  all  orders. 
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I  i  — —  a 

of  K  i  indeed,  this  sum  is  not  bounded  some  for  linear 

pi-oblems  of  this  tyi?e .  But  when  this  sum  1g  unbounded,  it  can  in 
fact  be  shown  that  the  general  boundary  value  problem  for  the  differen¬ 
tial  equation  concerned  has  no  solution,  and  that  therefore  in  this 
case  non-vanishing  solutions  of  the  associated  homogeneous  problem, 
l.e.,  elgen-ibnctions,  exist. 


5*  The  Boundary  Value  Problem  of  A  A ”0  .  In  order  to  show 
that  the  method  can  also  be  applied  in  the  case  of  differential  equa¬ 
tions  of  higher  order,'  we  consider  in  the  following  the  boundary  value 


problem  or  the  differential  equation 


£!<<  *  2 

lx*  ;x*V  V 


s  0  , 


We  seek  a  solution  of  this  partial  differential  equation  in 
our  region  (S  ,  for  which  the  function  values  and  their  first  deriva¬ 
tives  are  prescribed  on  the  boundary} . these  values  are  themselves 
defined  on  the  boundary  by  a  previously  given  function  (X  , y  ). 

In  this  case  we  assume  as  before  (p.  S^)  that  •f(x,y  )  le  continuous 
with  first  and  second  derivatives  in  the  plane  region  which  contains 
the  region  <3  . 

We  replace  our  differential  equation  with  the  problem  of 
solving  the ,  difference  equation  AA1X«0  for  the  grid  region  <3  , 


««  See  Courant-Hllbert,  Itethoden  der  matfaematischen  Physlk,  1, 
Chap.  Ill,  vbere  the  theory  of  integral  equations  is  discusseS 
with  the  aid  of  a  corresponding  alternative  principle.  See  also  the 
Q^ttingen  dlesertatlon  of  M.  v.  Koppenfels  (to  appear). 
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wherc  the  t'unction  'M,  assumes  the  tsame  values  as  the  previously  given 
(  X  j  y  )  at  the  points  of  the  boundary  line  r||  .  From  f  2  we 

know  that  the  boundary  value  problem  for  G  Jj  is  solvable  in  one  and  only 
one  way.  With  the  refinement  of  the  mesh  width  K  we  shall  show  that  in 
each  interior  sub-region  of  (p  this  solution  converges  toward  the  solution 
of  our  differential  equation,  with  all  difference  quotients  approaching 
their  respective  differential  quotients. 

For  this  purpose  we  first  notice  that  for  the  solution  of 

our  dif fei'ence  problem  the  sum 

«h' 

is  bounded  with  decreasing  mesh  width.  For  according  to  the  minimuin 
property  of  the  solution  of  our  difference  problem  (see  p.  12),  this  sum 
is  never  larger  than  the  corresponding  sum 

and  this  is  convergent  with  the  refinement  of  the  mesh  wiot;h  toward  the 


integral 


which  exists  according  to  our  assumptions. 


p'rom  the  boundedness  of  the  sum 


there  Immediately  follows  the  boundedness  of  jl  (A«)  , .  as  wed 

as  of  fv*  21  SI  )  and  In*  221  'Vi*  .  Bbr  arbitral^ 


Bbr  arbitral^ 


F-HMSRB  ?.'W'5PS?F 
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\/{f  there  exists  the  iheqaality 

k^'E.'Z  ur‘«cil*££  (15) 

G{^  G^^  *^k 

(see  (15)  p.  34),  If  one  substitutes  the  first  difference  quotients  of 
Mr  into  this  inequality  for  the  function  MT,  and  makes  use  of  the  sub¬ 
regions  of  Gjj  which  are  defined  with  these  difference  quotients,  there 
results  the  additional  Ineqjiality 

k*ZS  ^  e  k*  ZZ  )  +  “k  2:  (“/♦< ). 

6|i  <3fc 

where  the  constant  C  again  does  not  depend  upon  the  function  or  the  mesh 
width.  We  now  apply  these  Inequalities  with  luT*  It and  recall  in  this 
case  the  boundedness  of  the  sums  over  Pjj+  on  the  right  hand  side; 

the  boundary  sums  converge  by  definition  toward  the  corresponding 
Integrals  formed  with  f  (K  »  ^  Itiei^efore  from  the  boundedneas  of 

li'ZZCw*,****,* +u,*) 

^iT 

there  follows  the  boundedness  of 

k*£EK*0  .ni  k'ZSttV 

Sk  S|, 

We  now  substitute  successively  for  Mr  the  esrpftessions  Ah,  » 

the  inequality 

k*2S  +  ch*2£(Ai-)‘  (10) 

(fle<5  p.  30 ),  where  6*  is  a  subregion  of  G  containing  Q**  in  its  interior, 
all  of  which  satisfy  the  equation  A  WsO  ,  It  then  follows  that  for  all 
interior  subregions  (ff*  of  Q  f  the  sums 


s* 
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1  .e . , 


2 Z  ( ^ ,  K*  S  E  ( AH  i , 
s*  s* 

«re  teunded,  together  with  the  sums 

0^  ^ 


k”ZE(AH)‘. 

Sk 

Which  have  previously  been  shown  to  be  bounded. 

Finally  we  substitute  successively  for  tw  the  functions it  , 
Wi  ^ XXX  '  ^  inequality  (10), 

k”2£CAw-)‘,  1,..,  k'SZU'*.,)*..., 

®(.  S|* 

which  remein  bounded  as  Just  p:roved.  We  then  recognise  that  for  all 
subregions  the  sums 

2  X  (^XXX  ^jtx^  )  >  ^  ^^xyx'*'^xyy 

<3h* 

also  remain  bounded, 


From  this  fact  we  conclude  as  on  pp.  JO  ff,  that  we  may  choose  a 
partial  seq^eBoe  from  our  set  of  grid  functions  for  which  all  difference 
quotients  converge  uniformly  in  each  interior  subregion  of  G  toward  a 
continuous  limiting  function  in  the  interior  of  6  ,  which  are  their 
respective  differential  quotients.  \ 

We  still  have  to  show  that  this  limiting  function,  Which  obviously 
sstisfiCB  the  differential  equation  In  addition  fulfills  the 

prescribed  boundary  eonditions.  In  this  ease  we  satisfy  ourselves  as 
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above  by  expresHing  the  boundary  conditions  in  the  form^5 

Sj.  Sf. 

That  the  limiting  function  satisfies  these  conditions  may  be  shown  by 
applying  the  same  procedure  used  previously  on  p.  3^  the  function 
a^d  its  first  difference  quotients. 

0»i  account  of  the  unique  determination  of  the  solution  in  our 
boundary  value  problem,  we  can  now  recognize  in  addition  that  the  Indi¬ 
cated  convergence  properties  are  possessed  not  only  by  the  selected 
partial  sequence  but  the  function  sequence  li  Itself. 


II.  THE  IKPEHaOLIC  CASE 


f  1.  The  Equation  of  the  Vibrating  String 
In  this  second  part  of  the  article  we  are  concerned  with  initial 
value  problems  of  hyperbolic  linear  differential  equations)  and  we  shall 
show  that  under  certain  assumptions  the  solutions  of  the  corresponding 
difference  equations  converge  toward  the  solutions  of  the  differential 
aquation  during  refinement  of  the  mesh  width  of  a  basic  grid. 

We  can  most  slnqily  demonstrate  the  present  relations  with  the 
obvious  example  of  the  wave  equation 


(1) 


In  this  case  wa  restrict  ourselves  to  the  initial  value  problem  In 


It  is  not  difficult  to  show  that  boundary  values  for  the  function 
and  Its  derivatives  are  sufficient.  See  the  corresponding  observations 
by  K.  Friedrichs,  loc .  cit. 
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which  the  values  of  the  function  iA.  and  its  derivatives  on  the  straight 
line  "t-O  are  given. 

To  obtain  the  corresponding  difference  equation  we  construct  a 
uniform  square  grid  of  mash  width  H  iJii  the  X  ,  "t  -plane.  We  substitute 
for  the  differential  eqiaatlon  (l)  the  difference  equation 

tt  -  os  o  , 

j  '  O 

with  the  notations  of  pp.  5-4.  if  we  choose  a  grid  point  F®  ,  the 
difference  equation  in  this  case  connects  the  value  of  the  function  1/1 
at  this  point  with  the  values  at  the  four  neighboring  points.  If  we 
again  denote  the  four  neighboring  values  by  the  four  indices  1,  2,  5>  ^ 
(see  fig,  5)>  the  difference  equation  takes  the  simple  form 

u,  •♦■Uj  -  -u*  *0  .  (2 

In  this  case  the  value  of  the  function  K  at  the  point  P  itself  does 
not  appear  in  the  eqviation. 

'i  ■ 

We  imagine  the  grid  to  be  divided  into  two  different  partial 
grids,  as  indicated  in  fig.  5  with  circles  and  crosses.  The  difference 
equation  then  connects  with  one  another  only  tha  values  of  the  function 


Fig.  5 


Fig.  6 
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in  each  of  the  partial  grids.  We  therefore  restrict  ourselves  to  one 
of  these  two  partial  grids .  As  initial  conditions  ve  have  here  to 
prescribe  the  values  of  the  function  ii  on  the  two  grid  rows  =  and 
■ta  k  .  Nejct  we  state  e3Q)licitly  the  solution  of  this  initial  value 
problem j  l.e.,  we  express  the  value  of  the  solution  at  some  point  S  in 
terms  of  the  specified  values  on  the  two  initial  itjys.  One  can  isiasdl" 
ately  see  that  the  value  at  a  point  on  the  row  t  =  ak  is  uniquely 
defined,  merely  through  the  three  values  connected  with  It  on  the  first 
two  rows.  The  value  at  a  point  on  the  fourth  row  is  uniqjuely  defined 
by  the  value  of  the  solution  at  three  fixed  points  of  the  second  and 
third  rows,  and  therefore  by  certain  values  on  the  first  two  rows.  In 
general,  a  certain  domain  of  dependence  on  the  first  two  rows  will  belong 
to  a  point  S  }  one  can  find  It  by  drawing  the  linew  )f+t  e  const,  and 
X-t  «  const,  through  the  point  S ,  until  they  meet  the  second  row  In 
the  points  ei  and  ^  (see  fig,  .6).  We  thus  call  the  triangle  3*  ^  the 
diftteralnstion  triangle,  because  within  it  all  -values  do  not  change 
as  soon  as  they  are  specified  on  the  first  two  rows.  We  call  the  lateral 
lines  of  the  triangle  determination  lines. 

Denoting  now  the  differences  of  M,  along  the  determination  lines 
by  It  and  ft,  or  more,  precisely 

■w,  =  < 

«»=  , 

the  difference  equation  takes  the  form 

✓  ^ 

u,  *  w* 


This  means  that  a].ong  a  determination  Xlne  the  differences  in  the  other 

determination  direction  are  constant,  and  therefore  equal  to  one  of  the 

predetermined  differences  between  two  points  of  the  first  two  rows .  On 

the  other  hand,  the  difference  W.5  ~  is  a  sum  over  the  differences 

. 

along  the  determination  line  >  so  that  we  obtain,  by  application 
of  the  previous  remark,  the  final  formula 


in  which  the  notatlou  la  easily  understood, 

Nov  we  let  the  mesh  width  K  tend  toward  zero,  whereby  the  pre¬ 
viously  mentioned  values  on  the  second  or  first  row  converge  uniformly 
toward  a  twice  continuously  differentiable  function  ■f  (  )<  )>  »nd  the 

ii 

difference  qjiotlents  tend  toward  a  continuous  differentiable 

function  |9  (  X  ).  In  this  case  the  right  side  of  (?)  obviously  transforms 
uniformly  to  the  ej^iresslon  x 


If  if  converges  toward  the  point  ( t >  X).  This  is  the  well  know  expres¬ 
sion  of  the  solution  of  the  wave  equation  (1)  with  tbs  initial  values 
U(X,0)»'f(X)  and  ^<t/at(X,0)a  9  (X).  It  Is 

therefore  shown  that  the  Bolutlons  of  our  difference  equation  problem 
converge  toward  the  solution  of  the  differential  equation  problem  with 
decreasing  mesh  width,  If  we  let  the  initial  values  converge  in  tbs  above 
prescribed  scanner. 


-Uli- 


I  2.  On  the  Influence  of  the  Selection  of  the  Grid 
The  Regions  of  Dependence  for  Difference  and  Differential  Equation 

The  following  thoughts  are  obvious  in  view  of  the  considerations 
of  §  1.  Just  as  only  a  certain  part  of  the  initial  values  ere  decisive 
for  the  solution  of  a  linear  hyperbolic  differential  equation  at  a  point 
St  naiaely  that  "region  of  dependence"  out  out  by  chaise  ter  Istlca  through 
S)  there  Is  likewise  in  the  solution  of  a  difference  equation  at  the 
point  S  a  certain  dependence  region  which  one  obtains  by  drawing  the 
determination  lines  from  point  S  ,  In  §  1  the  directions  of  the  determi¬ 
nation  lines  of  the  difference  equation  coincided  with  the  directions  of 
the  cberacterlstiCB  of  the  diffe^ntial  equation,  so  that  the  dependence 
regions  also  agree  in  the  limit.  This  fact,  however,  depended  easentially 
upon  the  orientation  of  the .grid  in  the  (K  I't  ) -plane  and  was  based  on 
the  fact  that  we  heid  chosen  a  square  grid.  We  now  take  as  a  basis  a 
general  rectangular  grid  whose  mesh  width  In  the  -direction  (tljae-aash) 
equals  K  and  in  the  X  -direction  (space-mesh)  la  equal  to  xH  ,  with 
constant  X  •  region  of  dependence  of  the  difference  equation  >' 

^  completely  in  the  interior  of 
the  region  of  dependence  of  the  differential  equation  9 
or  will  contain  the  latter  in  its  interior,  according  as  X  <1  or  X>l  » 
From  this  follows  a  remarkable  fact;  In  tte  caae  X  < i  as  the 
mesh  width  K  decreases  toward  zero,  the  solution  of  the  difference  equa- 
tlon  cannot  converge  In  general  toward  the  solution  of  the  differential 
equation.  For  if  in  the  wave  equation  (1),  one  changes  the  Initial  values 
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/  \ 

^ *  »  ''' 

'  \ 


Fig.  7 

of  the  BOlutlon  of  the  differential  equation  in  the  nelghhorhood  of  the 
end  polotB  OL  and  ^  of  the  region  of  dei>endence  (see  fig.  7),  foirsula 

(4)  shows  that  the  solution  itself  also  ahanges  at  the  point  ( X  ), 

\ 

But  for  the  solutions  of  tlie  difference  equation  at  the  point  S  the 
initial  values  at  the  points  (H  and  ^  are  irrelevant,  Toeoause  they  are 
outside  of  the  region  of  dependence  of  the  difference  equation.  —  We 
shall  prove  that  there  is  convergence  in  the  case  >  I  in  f  J.  In 
this  connection  see  fig,  46. 

On  the  other  hand,  consider  as  an  exanple  the  differential  equation 
-  -0  (5) 

In  the  spatial  coordinates  X ,  ^  end  the  tljse  coordinate  I'  ,  and  replace 
it  by  the  corresponding  difference  equations  in  a  rectangular  grid.  Xn 
contrast  to  the  case  of  only  two  independent  variables,  it  is  now  iapos- 
sible  to  choose  the  nesh  spacing  to  that  the  regions  of  dependence  of  tbe 
difference  and  differential  equations  coincide,  for  tbe  region  of  dependence 
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of  the  difference  equation  is  a  rectangle^  while  that  of  the  differential 
equation  is  a  circle.  Later  we  shall  select  (see  §  h)  the  mesh  spacing 
BO  that  the  region  of  Influence  of  the  difference  equation  contains  the 
region  of  influence  of  the  differential  equation  in  its  interior,  and  we 
shall  show  that  convergence  then  occtu's. 

In  general  the  principal  result  of  this  section  is  that  one  can 
choose  the  grid  for  each  linear  hyperbolic  hoaogeneous  differential  equa¬ 
tion  of  second  order  such  that  the  solution  of  the  difference  equation 
converges  toward  the  solution  of  the  differential  equation  as  the  oesh 
width  tends  toward  zero.  (In  this  connection  see  7,8.) 


$  The  Faasage  to  the  Limit  for  Arbitrary  Rectangular  Grids 
React  nov  again  consider  the  mve  eq^ation 

n* 

end  seleqt  a  rectangular  grid  whose  time  mesh  width  Is  (t  and  whose 
space  mesh  Is  •  The  corresponding  difference  eqpiatlon  is 


i 

7 

'j 


where  the  indices  denote  the  central  point  P#  and  the  corners  , 

of  an  "elenentary  rhombus"  (see  fig,  8).  According  to  the 


A 

/V 

•  /  .  '  • 

/  > 

/  N 


-47- 


eqiiatlon  L(tt)=  0  we  are  able  to  describe  the  value  of  the  function  VL 
at  a  point  3  In  terms  of  Its  values  on  that  portion  of  the  two  rows 
iisiO  and  t«  K  ,  (see  fig,  6,  p.  4l)  obtained  by  drawing  the  two 
"determination  linea"  from  the  point  S  to  the  aides  of  an  elementary 
rhombus.  We  assume  the  initial  values  are  prescrlYjed  such  that  they, 
and  the  first  difference  quotients  formed  from  them,  converge  uniformly 
with  decreasing  mesh  width  and  with  fixed  X  toward  the  continuous 
prescribed  functions  on  the  straight  line  t  « o  ,  It  is  possible  to 
write  an  explicit  solution  of  the  difference  eqjiation  In  terns  of  the 
initial  values  (corresponding  to  (j)  In  I  1);  but  it  Is  not  so  simpl® 
to  perform  immediately  the  limiting  transition  to  vanishing  mesh  width. 

We  therefore  follow  another  course,  which  will, still  enable  us  to  treat 
the  general  problem.^^ 

We  multiply  the  difference  expression  L(<i)  ^9) 

write  the  product  under  consideration  according  to  the  following 
identities s 

) («,  -*«* ft<,)  s  ,  (7) 

24  In  the  following  compare  K.  Friedrichs  and  H.  Lewy,  "{Jber  die 
EindGUtigkeit  u.b.w,".  Math,  Annalen,  9B,  (1928),  pp,  I92  ff.,  where 
the  analogous  transformations  for  integrals  are  used. 


4. 
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We  thus  obtain 

Tjow  vn  flijm  the  TivriAuct  (9)  over  all  eleisentarj''  rhoabi  o*  a  dctcr=- 
mlnation  triangle  .  On  the  right  oide  of  (9)  the  aqneres  of  the 

differences  in  two  adjacent  elementary  rhombl  always  appear  with 
different  algns.  They  drop  out  in  the  aummstion  when  both  elementary 
rhombl  belong  to  the  triangle  $ ^  J  therefore  only  one  aum  remains  on 
the  boundary  of  the  triangle .  In  this  way  we  obtain  the  relation 

^  s 

Rare  H.  and  it  denote  dlfferencee  along  determination  directions  os  in 
§  1,  while  it  s-^^ote,0  the  difference  of  the  fimctlon  values  at  two 
neighboring  points  connected  by  a  line  parallel  to  the  "t  -axis .  The  sums 
are  to  extend  over  all  boxmdary  lines  consisting  of  two  parallsl  rows^  so 
that  a].l  differences  K  ,  U  ,  it  occur  once  and  only  once . 

Therefore  the  right  side  of  (lO)  disappears  for  a  solution  of 
L(ic)  =  e  .  Tha  sum  which  occurs  over  the  initial  rows  I  and  II  remains 
bounded  if  we  allow  the  mesh  width  ft  (with  fixed  K  )  to  decrease  toward 
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zero ;  in  particular  it.  tranafonn.';  Int-o  an  integral  of  the  given  function 
on  the  initial  line.  Consequently  the  sums  over  Soi  and  in  (10) 

also  remain  bounded.  It  is  now  that  we  mu.st  require  Jf  I  ,  (see  page  HU-) 
BO  that  I  “  is  non-negative  and  the  boundedness  of  the  individual  suras 

-  ^  T 

follows  which  we  can  imagine  to  extend  over  arbitrary  determination  lines. 
From  this  we  can  derive  the  "equl-continuity"  of  a  sequence  of  the 
grid  functions  In  all  directions  of  the  plane  (see  first  part  of  § 
because  the  values  of  %L  are  bounded  on  the  initial  line,  there  follows 
the  exiRtenoe  of  a  partial  seqjjence  which  converges  ■anlforuily  toward  a 
limiting  function  'U(J<,t). 

In  addition  to  the  function  1*  ,  its  first  and  second  difference 
quotients  also  satisfy  the  difference  equation  .  The  initial 

values  of  these  difference  quotients  are  expressed  by  means  of  the  equa¬ 
tion  L(«)«0  by  the  first,  second  and  third  difference  quotients  of  It 
in  which  only  points  on  the  two  initial  rows : I  end  II  appear.  We  require 
that  they  tend  toward  continuous  limiting  functions,  l.e.,  that  the  given 
initial  values  tt  (  X,0  ),  ^0  )  are  continuous  and  can  be  differen¬ 

tiated  three  and  two  times  with  respect  to  X  ,  respectively, 

25  If  S\  and  5*,  are  two  ixjlnts  separated  by  the  distance  S  ,  and 
if  one  connects  them  with  a  path  formed  by  two  lines  5,  S  and  $o^  , 
the  first  of  which  is  parallel  to  one  determination  direction  and  the 
second  is  parallel  to  the  other,  then  there  results  the  relation 
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We  can  now  apply  the  above  convergence  properties  to  the  first  and 
second  difference  quotients  of  IX,  instead  of  to  tX  itself,  and  we  can 


therefore  choose  a  partial  sequence  such  that  these  difference  quotients 
tend  uniformly  toward  functions  which  must  then  be  the  first  or  second 
derivatives  of  the  limiting  function  H(X  ).  The,  limiting  function  IX 
consequently  satisfies  the  differential  equation  *"  > 

which  corresponds  to  the  difference  equation  L(u)^(7  J  it  thus  describes 
the  solution  o.f  the  initiol  value  problem.  Since  this  solution  is  uniquely 
determined,  each  partial  sequence  of  the  grid  function  and  therefore  the 
sequence  Itself  converges  toward  the  limiting  function'. 


§  4.  The  Wave  Equation  in  Three  Variables 
We  now  consider  the  wave  equation 

,  »*K  _  _  »V  ^  g 

and  extend  the  reinarke  made  in  |  2  in  connection  with  the  region  of 
dependence.  The  region  of  dependence  of  the  differential  equation  (11) 
is  the  circular  cone  with  its  axle  parallel  to  the  i  -direction  and 
vertex  angle  «(  ,  with  +  J,M  «  s  s  .  In  a  reotangiilar  grid  ve  accord- 

,'12 

ingly  apply  the  difference  equation 


(11) 


Through  this  equation  the  function  values  44  are  connected  with  each 
other  by  the  points  of  an  "elementary  octahedron".  It  allows  the  function 
value  at  a  point  f  to  be  uniquely  e^^pressad  by  the  function  values  at 
certain  points  of  the  two  initial  planes  t  and  ■ts  k  •  We  obtain 


! 
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for  each  point  JS»  a  detemlnation  pyramid  which  cuts  two  rhonibl  from  the 
two  base  planes  aa  a  domain  of  dependence. 

If  we  let  the  mesh  width  tend  toward  zero  with  retention  of  its 
proportions  we  can  then  expect  convergence  of  the  sequence  of  grid  func- 
;  tions  toward  the  solution  of  the  differential  eqiiatlon  only  when  the 
determination  pyramid  contains  the  determination  cone  of  the  differential 
equation  in  Hs  interior.  The  simplest  grid  with  this  property  will  be 
the  one  which  is  placed  such  theit  the  determination  pyramid  touches  the 
determination  cone  from  the  outside.  Our  differential  eqjiation  is 
chosen  such  that  this  occurs  for  a  cuoic  rectangular  grid. 

In  this  grid  the  difference  equation  (12)  assumes  the  following 
form  in  the  notations  of  figure  10; 

here,  moreover,  the  function  value  Kj  at  the  middle  point  P  no  longer 
enters.  The  values  of  the  solution  on  the  two  initial  planes  are  the 
values  of  a  continuous  function,  four-times  dlrxerentlable  with  respect 
to  X,  /  »  t  .  , 

For  the  proof  of  convergence  we  again  use  the  method  developed 
In  f  5,  forming  for  the  solution  of  our  difference  equation  the  triple 
sun 

it’ SSI  2  L(1*)  »  0 

k 

which  extends  over  all  octahtdral  elements  of  the  determination  pyramid 

m 

enanatlng  from  a  pr  mt  g  ,  We  rei.ognlze  on  the  basis  of  an  almost 
literal  Interpr?  ..ntion  of  t!ie  previous  conclusion  that  the  values  of  the 
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Plg.  10 

IXmetion  ti  »t  the  intei’lor  pointe  of  the  determination  pyramid  drop 
out,  and  that  only  surface  sums  remain  over  the  four  lateral  double 
surfaces  F  and  over  the  two  initial  planes  I,  II  of  the  pyramid. 

Denoting  by  it*  the  difference  between  function  values  at  two 
points  which  are  connected  by  em  edge  of  an  octahedral  element,  there 
results  the  formula 

2£  (u')*  -  22.  («')*  =  °  •  ('•‘‘i 

F  rl 

which  is  to  be  s\ut#aed  over  all  planes  containing  the  differences  ^  ,  so 
that  each  such  difference  appears  only  once,^^  Since  the  double  sum 

The  grid  is  chosen  such  that  the  difference  of  %/L  between  the  two 
planes  f  no  longer  occurs. 


remains  tairx-ndeci  over  the  two  initial  planes,  as  it  transforms  into  an 
integral  o'oer  Initial  values,  the  sum  therefore  also  remains  hounded 
over  the  te-etenmiination  planes"  F  ■ 

Inila  cad.  of  to  'St  itself,  let  us  turn  our  attention  again  to  the 
first,  setomid  «nd  third  difference  quotients  which  satisfy  the  difference 
equation  (1!!;  5)  end  their  initial  values;  the  initial  values  are  themselves 
expressedijc'yiaeons  of  (ij)  by  the  first  througti  fourth  difference  quo¬ 
tients  foniCTied  d’rom  the  valuts  on  the  first  two  initial  planes.  If 
Uraur^i)  oac  of  these  difference  quot,ient8  up  to  the  third  order,  we 
then  knovlrt  hat  the  sum  remalna  bounded  over  one 

determinatlo on  ]plane  F  .  By  exactly  the  eaine  consideration  which  ve 
applied  int^ths  .first  part  of  J  la,  it  follows  that  the  function  li  and 
its  firstarx-jd  *»oeond  difference  quoti^ents  are  equl- continuous.  Therefore 
there  exJitme-as  sequence  of  (nssh  widths  which  deereasee  toward  *ero  such 
that  thosaocxpxeas ions,  which  are  bounded  In  the  beginning,  converge 
toward  coitt-luoua  limiting  functions,  and  indeed  obviously  converge 

toward  tho » 80l%itlon  of  the  differential  aquation  Including  its  first  and 

11  '  ■  ,  ■ 

sedond  dfli!lv#'-vatd.vs8,. which  follows  eseaetly  as  In  |  J. 


-5l+- 

APPEWPIX 

SUPPIJMECTS  AND  GENERA-LIgATIOHS 
§  5*  ExflHnple  of  r  Differential  Equation  of  the  First  Order 

We  have  seen  In  §  9  that  under  soiie  circumatanees  the  domain  of 
dejnendenne  of  the  differential  equation  oonotltuteo  only  a  part  of  the 
domain  of  dependence  of  the  difference  equation,  and  therefore  the 
influence  of  the  remaining  region  drops  out  in  the  limit.  We  can  expli¬ 
citly  demonstrate  this  phenomenon  with  the  example  of  the  first  order 
differential  equation  ^0  if  we  substitute  for  it  the  difference 

equation 

Written  in  the  notetionn  of  fig,  5  (p,  111)  it  reads 

;  (16) 

f.  ■  ' 

This  dlffersnee  aq,uation  again  connects  only  the  points  of  a  partlajl 
grid  with  one  another,  The  initial  value  problem  consists  of  specifying 
the  function  tt  as  the  values  -f  which  a  continuous  function  ■f  (  X  ) 
assunes  at  the  points  Xa  2ik  on  the  row  ^  sO  . 

Consider  the  point  S  on  the  t  -axis  at  a  distance  2hH  •  It  Is 
easy  to  verify  the  representation  of  the  solution  1C  at  i?  as 


The  sum  on  the  right  side  tends  toward  the  value  with  rafinsmant  of 
the  mesh  width,  i.e,,  as  .  One  can  conclude  this  from  the  conti¬ 

nuity  of  4 (  X )  and  the  behavior  of  the  binomial  coefficients  with 
increasing  n  .  (See  the  following  paragraph.) 
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§  6.  The  Eqjiation  of  Heat  Conduction 
The  difference  equation  (l6  )  of  §  5  also  be  considered  an 
analogue  of  a  quite  different  equation,  namely  the  equation  of  heat 
conduction 

2 


=  0  .  (18) 


In  any  rectangular  grid  the  corresponding  difference  equation  reada 

where  I  Is  the  tlioe  aosh  and  tbo  space  mesh.  In  the  limit  of 
vanishing  mesh  width  the  dlfferenos  equation  laaintalns  its  form  only 
when  1,  decreases  In  proportion  to  It^  ,  In  particular,  If  we  eat 
I  ■  K*  ,f  the  value  Hq  then  drojm  out  of  the  equation  and  there  WBults\^ 
the  difference  eqpjiation 

■u,  =  tlLlliS  ,  (16) 

i 

whose  solution  has  been  given  by  formula  (IT): 

<*  oA )  ”  i:  { »*" )  Ki  ■  h7) 

A  point  t  of  the  X  -axis  is  alvays  denoted  vitb  deoreasing  mesh 

/  //' 

width  by  the  index 

21  a  i  •  '  ,  (SO) 

9, 

The  mesh  width  is  related  to  the  ordinate  “t  of  the  source  point  by  the 
equation 

Zn\^  s  i  .  (21) 

We  shall  now  examine  what  re  suits  from  fomala  ( 17 )  when  fv  tends 
toward  aero,  l.e.,  ae  Ifi  tends  tOTrard  Infinity.  By  application  of 


(22) 


forniulft  (21)  we  may  wi’ite  equation  (l?)  In  the  form 

”  /  2« 


it  i  ) 


U-H 


Z‘Z 


an 


For  the  coefficients  of  we  use  the  abbrevi¬ 


ation 


Or 


*Ft  ”  )  ■ 

Here  wa  shall  calculate  the  limit  of  this  coefficient,  which  one  usually 

i: 

determines  with  the  aid  of  Stirl))iig*D  formula,  by  interpreting  the 

y' 

function  9  jn,  (  I  )  ®“  *  ordinary  diff«j‘enee  eqtuation, 

and  proceeding  to  the  limit  for  vanishing  nesh  width  k  and  thus  to  the 
differential  equation.  One  finds 

as  the  difference  eqiiatlon  by  writing  )  Instead  of  $ )• 

^(9k««M-JiC»»)  -  ■■3k<S)d^b;i  • 

In  addition  | )  satisfies  the  nonsAlisation  condition 

^  0i,(<)-ah  a  2 )rr . 

This  BuiB  is  to  eictend  ovai*  the  dotsaln  of  dependsnoe  of  tbs  difference 
equation,  which  in  the  limit  as  k-^  0  occupies  the  whole  X  -seels. 

One  can  see  by  simple  considerations  that  9kOi)  converges 
uniformly  toward  the  solution  3  ( ^  differential  equation 

gVx)  =?  “  r 

"ti 

with  the  auxiliary  condition 


•o 

I  $Cx)<(x  ==  zft  . 


-•-0# 
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By  proceeding  to  the  limit  there  then  arleea  from  fonnula  (22) 

the  well-known  solution  of  the  equation  of  heat  ctmduetibn. 

The  cons iderat ions  of  this  paragraph  are  directly  applicable  to 
the  case  of  the  differential  equation 

wau  ^  aV  ^ 

etc.  in  Bore  Independent  variables. 


§  7,  The  General  Linear  Homogenous  Differential  Equation 
of  Second  Order  in  the  Plane 

V.  ,  '  . 

Let  UB  consider  the  differential  equation 


(2» 


The  coefficients  are  twice  continuously  differentiable  with  respect  to 
)C ,  “t  ,  while  the  initial  values  pn  ths  straight  line  teO  are,  three 
tlmss  oontlnubusly  differentiable  with  respect  to  K  .  Ve  replace  the 
differential  equation  by  tha  difference  equation 

LCo)*  +  ■'■i'*'"®  W 

in  a  grid  with  the  time  meeh  width  H  and  the  space  mesh  width  xk , 
such  that  in  the  nel^borbood  of  that  part  of  the  initial  lltte  under 
consideration  I  >€  >0  holds  for  our  constant  X  i  we  choose  the 
Initial  values  as  in  f  J,  (See  p,  1*9.) 

Iibr  the  convergence  proof  ve  again  transform  the  sum 
L*  tr  c  «  ill  “  I  / .  V 


2  L(«) 


by  application  of  the  identities  {7)>  (6).  In  addition  to  the  sum  (see 
(lO))  over  the  boundary  of  the  triangle  (see  fig.  6)  there  nov 

appears  a  sum  over  the  vhole  triangle  Se(^  ,  whose  absolute  value  can 
be  estimated  with  the  help  of  the  Schwara  inequality  as 

where  the  constant  C  ^*8  not  depend  on  the  function  U,  on  the  mesh  ^ 
width  k  ,  or  on  the  point  -6  in  #  nertaln  neighborhood  of  the  initial 
line. 

Here  we  can  almilarly  estimate  the  value  of  h  'It*  as 

C.I.‘2£(|)%C.kZv‘. 

/H  t*  .  . 

Where  what  has  been  said, for  C  vlll  be  valid  for  the  constants 

He  thus  obtain  an  Inequality  of  the  fom 

,  t » .  ») 

wt(ere  D  is  a  fixed  constant  for  all  the  sums  over  the  initial  etraigtit 
line^  for  all  points  S  and  oesh  widths  K  • 

Starting  from  the  initial  line  we  nov  choose  as  the  point  S  of 
our  triangle  the  points  S0  t  S^p  •••»  ^  in  turn  on  a  line 

parallel  to  the 't  -axis.  By  suanatlon  of  the  inequalities  corresponding 


For  the  proof,  see  the  related  Inequality  at  the  bottom  of  p.  49* 
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to  (25)  we  obtain  the  ineqviallty 


ss^<f( 

•''■z§[K'-t‘)sr*gRn 

✓  \ 

Reoftlling  that  one  can  ez^iresa  the  difl^rence  14  or  14  la  tanw  of  thii 
two  dlfferesaoB  it  end  a  difference  14  or  It ,  reBpectively^  It  foXXovs 
that  we  can  then  reduce  the  left  side  of  (a6)  if  we  aubatitute  for  it 

We  now  reatriet  ouraeXvea  to  •  naiflhborhood  t  ivk  of  the  inltiaX  Xine 


in  which 


C^  -hK  C^  «  Cg  > 


BO  va  obtain  fros  (s6) 


c,V 


Fro«  thla  aBSlifeaBion  of  the  bwadidaeii  of  the  Xaft  aide  of  (27),  it 
foXlowB  froB  (85)  that  ;  . 


ifl  a|.ao  bounded,  ftoM  which  thara  fOXloWB  the  aqul-contlnulty  of  14  a$ 
in  J  3, 

Inatead  of  applying  the  ine(juallty  (25)  to  the  function  lo  itealf 
we  apply  it  to  the  first  and  second  difference  quotients  l^r ,  which  also 
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aatlsfy  difference  equations  whose  second-order  terms  are  as  in  (24). 

In  the  additional  terms  derivatives  of  ii  which  cannot  be  expressed  by 
can  indeed  still  occur,  but  the  sum  of  their  squares  over  a  rectangular 
area  nailtiplied  by  ,  ean  already  be  assumed  to  be  bounded.  But  la 
this  qase  ve  may  apply  to  this  difference  equation  for  Kr  t^jie  same  con¬ 
clusions  which  have  been  applied  earlier  to  <U  .  Therefore  ve  can  infer 
the  equi -continuity  and  boundedness  of  the  functions  1i  and  their  first 
and  second  derivatives,  which  therefore  possess  a  partial  sequence  which 
converges  uniformly  toward  the  soli^\tlon  of  the  Initial  value  problen  of 
the  differential  equation.  Trm  the  uniqueness  It  again  follows  that  the 
function  salience  itself  convorgea.  n  . 

I  In  this  case  of  course  we  Bust  assuaa  that  the  dif ferenoe  dHotieots 
IIP  to  the  third  order  oonvsrgs  unlfqsaly  toward  oontlnuous  llaitlBg 
fUnotlona®®  on  and  hetvien  the  initial  rows. 


I  8.  The  Initial  Vhlue  Problea  of  ah  ^arbitrary  aypet^llo 
Uhtar  Differential's^ 

VTa  shaii  now  stow  that  the  Mthodi  davaloped  above  ar»  adequate 
to  solve  ths  Initial  value  pfoblea  of  ib^h^itnry  lioelur  bc^ 
hypex^Iio  iiUffirential  aquation  of  s  orisr,  Xu  this  cast  it  is 
eufflclent  to  rsstriot  ourselves  to  the  case  of  three  variables.  The 

'  v  '  „■  "  ■ 

train  of  thought  aay  be  lanedlately  applied  to  several  variablas.  One 


can  easily  see  that  the  most^'^gsneral  problem  of  this  kind  can  be  reduoed 

■  A  '■  ...  ■  - 


_  V  , 

^  This  assuBiption,  ^end  those  coaoemlug  the  differentiability  of 
the  coefficients  of  the  differential  equation  and  the  boundedness  In  a 
sufficiently  sas.ll  neighborhood  of  the  initial  lines,  can  be  relaxed  in 
particular  coses. 


-6l- 


by  a  trans format ion  of  variables  to  the  following:  to  find  a  function 
which  aatlsfiee  the  differential  equation 

'‘♦it  ""  (23) 

and  yixiQhf  with  ita  ilrat  darlvatlvaa,  aaauaaa  presorlbed  values  on  the 
plane  taO  ,  in  this  ease  the  coefflelente  may  be  functions  of  the 
variables  ^  ^  ,  t>  should  aatiafy  the  coaditioas 

A.  >  o  ,  c  >0  •  Ac-b*  >  o  . 

Ve  also  aasuae  4hat  the  aoefflolents  are  three  tinea  continuously 
'differentiable  with  respect  to  X  i  t  initial  values 

Y<are  four  tlane  differentiable  and  14 «  three  tunes  differentiable  with 

“  I'l  ■  .  ...... 

\iieiipict  to  X  f  ^  .. 

tf#  BOW  aaauae  the  eoordinatea  X  and  ^  are  rotated  around  a  point 
is  the  ihiibial  plaiw  suoh  that  b«0  .  Than  in  a  certain  neighborhood  0 
of  this  point  thav  eondltlona 

'  \4-|b|>«^»  c-jbvxp 

m  itiraly  fttlfiUeidt  tfe  liait  oiir  oonaidaratioiiie  to  thia  mlgbbdrheod.  /, 
Vtf  thra  eln  oh^e  i  th^  tlMi  oontinuoualy  dlffereniiabit  funotlon 
(1>P- •uoh-tbat- 

,  a  -  4.,  V  i-  , 

.  -  ,1  ■  4  -,lb|  ■;  ■ 

ig  Tslid  with  constant  $  .  We  aan  then  put  the  differential  eiiaation  in 
the  form 

'll 

-  (o.- 4 )  *  (c- 4  J j  (H+ b)(‘U^x  +  + 'Wjij, ) 

■  \ 

“‘^(4-b)C;Uxy»t'K^4H^9)  + +8u  =  o.  (50) 


/ 


(89) 


/;  . 


■  ,,  r 

-  ,7 
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We  now  construct  the  grid  r>f  points 

is  ill ,  x+y  a  tnxK  ,  x-y  *  hick  (  l,m,h  a  o, »,  a, " * ) 

in  the  region  and  suhetltute  for  equation  (JO)  a  diff^l-ence  equation 
L(  <4)aO  in  thlB  grid.  For  this -puiTpbae  ve  adjoin  to  each  grid  point 
the  following  adjacent  points;  ^[5ie  ^Ints  and  ,  which  oob» 
fioa  di^aplacenenta  k  and -k  »  respectively,  in  the  dlieetion  > 

pf  the  i  -axis}  ,,aod  the  points  Pj  ,  ...>  P^  ,  which  lie  in  the  saoe  plane 
as  P^  parallel  to  the  (X  ,  ^  )-plane;  see  fig.  11.  These  points  fora  an 
^^octahtf^ral  aleaent'* ..with  the  oonwsr  points  ,  P^  ,  Pj  >  P*,  t  ^3  i. 

t."  •  ■ 

^  --X-  - 

fov  each  grid  point  r*  vhish  il ’'inside  6  ve  replace  the  stoohd 
difiterentlal  q^tients  comring  ih  (^)  hy  difference  quotients 


/7e 


♦I 


*0 


9 

3 


•3 


pig."ii''; 

about  P^  in  tbs  octahedital  eleaent  in  the ;  fOlloving  way. 


We  replace 


(*.  h 

W  jijj.  (■«. , 


t<K){ '"  *  '•*  ■*■  *^Sh> 


W«  nil)«tit:i;t«  oorreBpondlsg  difftr«ae«  (motieati  is  o«t*lMldr»l 
•iMBiftt  t<at  the  fltst  dlfibnatil*!  qaotlentB  oooweviag  in  (30) »  W«  igl'ra 

.  •■  ,  'i  ■  „  -  " 

ttM  flonffioicnti  in  ih«  diffttWDO*  f^tinn  the  valuBia  vhlQh  «r«  taken 
'^t>y  the  ooafflolantB  of  tha  ^ifftrintial  aqaaUojii  at  tiw  point  r*  • 

V t  ie*  «ttd 
IVmotlonrvaluai  pwiarilJtd  ad  ttait  tbay  -aonvajp^  toward- tba‘pifiimerik<B 
"initial' ^uaa  at  .iJ  ad  wits  .si'ftolia®t!lofitiM»::«aah'_witoh.rwhliiL::^ 

;  tht.  apa8i-''to  tMrtlwM^';v^ 

■  '  '■(('  ■  lU'.'/i'  '1  _ 

diffawnce  ^otianta  up  to  tha  fbayrtS  PJ^r  feri^ 

tso  and  ta  K  ahould  mif^raiy  oojiy^  toward  tne  oprnnippndlag 

praaerlbed  diffarant laX V quotlanta t 

■■■'’"■  •  ,  ,.  '"  ■■  .  '  I 

Tha  aolutlon  of  the  dlffareacaiaqjpiation  LCit)*®  at  a  point  la 

■.  ■•''  ■  ■  if 

unlQuaXy  datamlnad  by  tha  valuaa  on  the  two’botton  plauaa  of  tha  datar^* 
nlnation  pyraald  paaaing  thr^gh  it.  \ 


t 
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For  the  proof  of  convergence  we  form  over  all  elementary  ootahedrona 
of  a  determination  pyramid  the  sum 

;  k’ s  ui.) 

and  transform  it  hy  means  of  the  identities  (7)/ (8).  In  this  way  there 

'  "I »  '  " 

arises  one  spatial  sum  itoltlplied  toy  It  ,  which  la  «iuadratlc  in  the  first 
difference  quotients,  and  also  a  siw  multiplied  toy  \  Jj*  over  the  lateral 
^tole  planes,  in  vhioh  appear  the  sqparea  of  aU  difference  quotients  of 
tha  type  occurring  on  and ‘between  the 

double  planes}  their  coeffieiente  are  greatar  than  a  find  pbiitlve 


ooheti&t  beeauie  of.  if  chooie  moreover  e  su:IH'loi«ntly  eaeii  ratio 


I 


»• 


I  ■  ((,'  ■ 

X  Of  tine-  to  epace-oeeh  vidths , 

prow  here  we  can  proceed  in  tl»  line  iiaehion  .ill  117/  and 
can  prove  that  the  solution  of  our  difference  eqtmtidtt  convargaa  toward 
the 


/!' 


t// 
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